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Abstract 

In this paper we present a regression based model for day-ahead electricity spot prices. We estimate the 
considered linear regression model by the lasso estimation method. The lasso approach allows for many possible 
parameters in the model, but also shrinks and sparsifies the parameters automatically to avoid overfitting. Thus, 
it is able to capture the autoregressive intraday dependency structure of the electricity price well. We discuss in 
detail the estimation results which provide insights to the intraday behavior of electricity prices. We perform an 
out-of-sample forecasting study for several European electricity markets. The results illustrate well that the efficient 
lasso based estimation technique can exhibit advantages from two popular model approaches. 


1 Introduction and motivation 

In day-ahead electricity spot price forecasting there exists a vast range of differently structured models. A recent 
overview with a classification is given in Weron (2014). 

Most of the day-ahead electricity spot prices are traded on an hourly (or half-hourly) basis. This leads to the 
obvious fact, that we have 24 (or 48) observations a day if we ignore the clock change issue. Subsequently we only 
consider hourly traded prices, but for half-hourly data the adjusted facts hold as well. All electricity price models 
rely heavily on this fact of having 24 observations a day. These models can be divided into two classes independent of 
whether the model is designed e.g. for point or for probabilistic forecasts. 

Modeling approaches which belong to the first class construct 24 separate models, one model for each hour of the 
day. These models can be e.g. a standard time series method (e.g. Paraschiv et al. (2014), Nowotarski and Weron 
(2014),Maciejowska and Weron (2015)) or a regime switching model (e.g. Karakatsani and Bunn (2008b)). In general, 
the model approach is usually motivated by the fact that there is one auction each day (mostly in the morning or at 
lunch time) for the next 24 hours^. Hence, we observe the 24 traded prices at once. Thus, modeling a 24-dimensional 
process by 24 small models is an obvious way to proceed. 

The second class contains electricity spot price models that propose a single model for the full observed data. Here 
a 24 hours ahead forecast is computed. This approach is motivated by seeing the price process as one time series that 
evolves over time as related processes like electricity load or renewable energy generation. As in the first model class it 
can be of various model types, like time series based approaches (e.g. Haldrup et al. (2010), Cruz et al. (2011), Keles 
et al. (2012), Paraschiv et al. (2015) or Ziel et al. (2015a)) or machine learning models (e.g. Conejo et al. (2005), 
Singhal and Swamp (2011) or Chaabane (2014)). 

From the applied point of view both model approaches are valid as they can have a good forecasting performance 
in particular situations. Nevertheless, the first one seems to be more plausible, as it matches more directly the data 
generating process. However, in the first class with 24 models for each hour we have usually less observations and 
consequently 24 simple and low parametrized models (see e.g. Weron and Misiorek (2008) or Maciejowska et al. 
(2015)). More crucially, most methods assume the same model structure for all 24 hours. The model structure is 
often assumed in such a way, that the price at hour i depends on the price of the previous day at hour j, but not on 
a different hour J. Many models do not allow for such a dependency to simplify the model structure, even though 
there might be strong dependency present. Here models of the second class can have an advantage. As they use more 
observations for the estimation the resulting model tends to have more involved parameters that can theoretically 
capture the dependence, if the model design allows for it. Interestingly the vast majority of machine learning models 
(esp. artificial neural networks) fall into this second model class that estimates one big model, see e.g. Conejo et al. 
(2005), Singhal and Swarup (2011), Catalan et al. (2011), Shafie-Khah et al. (2011), Hu and Taylor (2014) or Abedinia 
et al. (2015)). It seems that these estimation techniques rely much more on the fact of having a lot of observations 
than time series based approaches. 

To better motivate and understand this dependency problem we introduce Pd,h as electricity price at day d and 
hour h. Throughout the paper we ignore we clock change issue, and interpolate the missing 2am hour in March and 
average the double occurring 2am hour in October so that there are 24 observable prices each day. 

Now we have a closer look at the linear dependency structure between the price Pd^h at day d and hour h and 
the price Pd-i,i at the previous day d — 1 and hour 1. The corresponding sample correlations of the correlations 

^Once a year for the last Sunday March only 23 and for last Sunday in October 25 hours due to daylight saving time 
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CoT{Pd,h, Pd-i,i) are exemplary visualized in Figure (1) for two selected markets, namely the EPEX spot price for 
Germany and Austria, and the APX spot price for Netherlands. 
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(a) EPEX spot price for Germany and Austria 
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(b) APX spot price for Netherlands 
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Figure 1: Sample correlations of Cor{Pd^h, Pd-i,i) for h, Z = 0,..., 23 and two markets from 17.12.2009 to 12.08.2014. 


The highlighted diagonal elements measure the linear dependency between the electricity price at hour h to its 
previous day value at the same hour h. We observe relatively high correlations on the diagonal. But, only for some 
hours in the night, the late afternoon and the evening they are the largest of each row. For most of the hours from 
midnight to the afternoon this is not the case. They exhibit their largest correlation with hours in the late evening. 
The most distinct lagged cross-hour dependency concerns Td,o- It depends much more on the last hour price Yd- 1^23 
of the previous day than on the same hour price Yd-ifl. 

It is clear that this complex dependency structure should be covered by an appropriate spot price model. As 
mentioned before, most of the models of the first class do not capture this dependency structure, whereas the second 
one in principle does. But, the second class that uses one big model has the problem, that e.g. the forecast Fd+ 1,23 
will depend on the previous hour predictions Yd+i ^,... ,Yd+i^ 22 - Thus the prediction error ^d+i,h = Yd+i^h — 
tend to increase in h as we have to compute 24 hours ahead forecast. And even though models of the second class can 
capture the cross-hour dependency, we observe that the full dependency pattern of Figure 1 is time-varying. Thus, 
it is required for models of the second class to utilize time-varying parameters for the autoregressive effect. This is 
usually not included in many models of this model class, see e.g. Conejo et al. (2005), Singhal and Swarup (2011), 
Catalao et al. (2011), Shafie-Khah et al. (2011), Hu and Taylor (2014) or Abedinia et al. (2015). 

Interestingly, when we compare forecasting performance of relatively simple models of both model classes, models 
of the second class perform better for the first half of the day, whereas models of the first class are better in the second 
half of the day. This observation is also part of the results of the empirical section. Nevertheless, it is not clear which 
modeling method is superior in terms of forecasting performance. 

In this paper we introduce a simple linear model approach that falls into the first model class. It contains 24 simple 
structured models and captures the mentioned dependency pattern directly. Thus, it joins the advantages of both model 
classes which leads to an improved overall forecasting performance. The model captures the full linear autoregressive 
dependency structure as seen in Figure 1, especially time-varying cross-hour dependencies. Furthermore, it also allows 
for a time-varying day-of-the-week effect that will be discussed more detailed in the next section. In electricity price 
forecasting literature such flexibility is a clear novelty even though there are models available that cover the same 
feature. For example, an autoregressive principal component analysis based approach as used in Maciejowska and 
Weron (2013) and Maciejowska and Weron (2015) can cover this behavior. However, there the impact is not modeled 
directly but indirectly through the underlying dimension reduction approach. 

We restrict the model to these two modeling aspects (time-varying cross-hour dependencies and time-varying day- 
of-the-week effects) to highlight the advantage of the modeling and estimation approach. Especially, we do not focus 
on well known features like volatility clustering, non-linearity effects (esp. price spikes), long-term trends, structural 
breaks, public holidays effects, effects of other external regressors like renewable energy feed-in, or the day-light saving 
time effect. However extensions to all these directions are possible. We discuss this issue briefly in the model section 
and final discussion. 

The crucial part of the estimation relies on high-dimensional statistics. We utilize the lasso (least absolute shrinkage 
and selection operator) method introduced by Tibshirani (1996), as a regularized estimation technique. Due to the 
shrinkage properties of the estimators we can handle models with many potential parameters. As the estimation 
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procedure has the sparsity property, many of the potential parameters are not included in the finally estimated model. 
This property is great to avoid the overfitting problem. 

The considered estimation technique is very fast. The estimation and forecasting of the 24 models is feasible in 
a few seconds without big effort. Thus, the model is also suitable as being an element of a forecast combination 
approach, see e.g. Bordignon et al. (2013) or Maciejowska et al. (2015). For such applications fast estimation methods 
are clearly preferable. 

The lasso based approach is not completely new in electricity price forecasting literature. It is used by Ziel et al. 
(2015a) and Ludwig et al. (2015) for electricity price forecasting. Both models use the variable selection property, but 
consider a model of the second class and share the mentioned disadvantages. The approach of Ziel et al. (2015a) uses 
time-varying coefficients to capture the mentioned time-varying changes in the dependency structure. However, the 
time-varying coefficients are modeled using B-splines, that only allow for smooth changes over the day. In contrast, our 
model belongs to the first model class and allows for more flexibility and abrupt changes in the cross-hour dependency. 
Additionally, Ziel et al. (2015a) incorporate the impact of load and wind and solar net feed-in such as volatility 
clustering, public holiday and daylight saving time effects. In Ludwig et al. (2015) the parameters are assumed to 
be constant over time, so the lasso estimation works just as features selection and ignores time-varying dependencies. 
They use wind speed and temperature data to model exogenous impacts. In recent research of Dudek (2016) the lasso 
approach is applied to short-term electricity load forecasting. 

The paper is structured as follows: In the second section we introduce the time series model and estimation method. 
In the third one we discuss several benchmarks for the empirical study which is performed in the fourth section. For 
the out-of-sample sample forecasting study we consider 10 different hourly European day-ahead electricity spot prices. 
We close with a summary and some conclusions. 


2 Time series model and estimation technique 

The time series model that we consider is an autoregressive regression approach, that is similar to the one used in 
Weron and Misiorek (2008) or Maciejowska et al. (2015). In Weron and Misiorek (2008) the electricity price at hour h 
depends linearly on the price at the same hour lagged by 1, 2 and 7 days, such as dummies on the weekdays Sunday, 
Monday and Saturday. However, the choice of the lags 1, 2 and 7 such as the selection of the weekday dummies is the 
same for all 24 hours and it is very restrictive. Concerning the weekday structure it is advisable to take a closer look 
at the weekly mean price of the electricity price. The weekly sample mean of the German/Austrian EPEX spot price 
and the Dutch APX spot price is given in Figure 2. 



Hour 


Hour 


(a) EPEX Germany and Austria 


(b) APX Netherlands 


Figure 2: Weekly sample mean of Pd,h for two markets. 


There we see the 7 x 24 = 168 sample mean values grouped by the 7 weekdays and 24 hours of the day. We can 
observe that the sample mean of the Sunday and Saturday values is clearly different from the mean on working days. 
On the working days from Monday to Friday the mean behavior is relatively similar. But for the morning hours on 
Monday and the afternoon and evening hours on Friday we observe clear deviations. This transition towards and 
from the weekend is known for the German/Austrian EPEX and modeled by Ziel et al. (2015a) by an approach of 
the second model class using time-varying parameters. Most importantly it suggests that the impact of the weekday 
















2 Time series model and estimation technique 


4 


dummies in an autoregressive regression model might be different for different hours of the day. Hence, we require 
flexibility in the model that is able tackle these possible changes in the day-of-the-week dependency structure. 

For the time series model that we introduce for the electricity price Pd^h we assume that the processes {Pd,h)dez si's 
stationary for each hour h. Furthermore, we do not directly model Pd^h itself, but its zero mean version Yd^h = Pd,h — ldh 
with = ]E(Fd,/i)- We can estimate fdh easily by the sample mean. 

To model the mentioned weekday effects we define the weekday indicators 


Wk{d) 


1 ,yv{d)=k 

0 ,W{d)y^k' 


Here >V(d) is a function that gives a number which represent each day of the week. We use without loss of generality 
the enumeration A: = 0 for a Sunday, fc = 1 for a Monday, ..., fc = 6 for a Saturday. As we want to point out the 
effect of the weekdays, we propose two models, one without the weekday effects, the other one with possible weekday 
effects. 

For each hour h the two considered time series models for Yd^h are given by 


24 

Yd,h = ^ ^ (t>h,k,iYd-k,i + £d,h and 

l = l kGX{h,l) 

24 6 

Yd,h = E E 4>h,k,lYd-k,l + ^ i’h,kWk{d) + ed,h 

l = l kGX{h,l) fc=0 


( 1 ) 

( 2 ) 


with parameters (j)h,k,i and V'/i./o 21(/i, 1) as an index set of possible lags and £d,h as error term. We assume that the error 
process {£d,h)dei. is uncorrelated with constant variance cr^. The parameters 4’h,k,i model the linear autoregressive 
impact of the past prices and iph,k a weekday effect. 

The index sets T{h, I) are very important for the model. They specify the possible model structure. Of course larger 
sets I{h, I) increase the probability of overfitting, even though the used regularization technique reduces this problem 
a lot. If we choose I{h, I) too small we might miss important information in the dependency structure. Therefore we 
should choose T{h,l) reasonably, in the sense that we should only include those lags where we, as a modeler, think 
that a dependence is plausible. In our situation we choose 


T{h,l) 


{!,..., 36} ,h = l 
{!,..., 8} ,hj^l 


(3) 


Thus, we include a possible relationship between Yd^h and the previous 36 days (five weeks plus a day) at the same 
hour h, but also a dependence to the past 8 days of the electricity prices of other hours. 

Furthermore, model equation (2) (also (1)) is a regression model, so it is very modular. It is extremely easy to 
add new regressors, like lagged electricity load or renewable energy feed-in, seasonal cycles, trends, or event dummies. 
These additional regressors can be used to cover some of the ignored stylized facts, like long-term trends, structural 
breaks, public holidays and day-light saving time effects. Even non-linear impacts in Yd-k,h like in a threshold-AR 
process are possible, as long as the relationship is linear in the parameters. Another direction of extensions concerns 
the volatility of the residuals. Ziel et al. (2015a) and Ziel (2015) show that an extension to AR-ARCH type processes 
is possible as well, using an iterative lasso estimation technique. 

For the estimation of the parameters in (1) or (2) we will use the corresponding ordinary least square (OLS) 
representation. We denote it by 


Yd,h — X’d^hfdh + ^d,h- (4) 

with the vector of regressors Xd,h = {Xd,h,i, ■ • ■, Xd,h,m,,y of length mh and parameter vector (3^ = , I3h,mhy 

of length mh- Furthermore, we require for the lasso estimation technique the standardized version of (4). Thus, we 
introduce it with 

Yd,h = ^d.hf^h 4" ^d,h- (5) 

Here and the elements of Xd^h are scaled in such a way that they have a variance of 1. In the empirical part we 
scale using the corresponding sample standard deviations. If we have we can calculate the parameter vector /3^ of 
(4) easily by rescaling. 

The mentioned lasso estimation technique is a penalized regression approach that penalizes the parameters on the 
sum of the absolute values. Let D be the number of observable days. Then the lasso estimator f3h for f3h is given by 

^ D _^ 2 

3?, = argmin V -HA?,V|/3i| 

d=i ^ 


(6) 
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where A/i > 0 is a tuning parameter. Note that for = 0 we receive the standard OLS estimator. For extremely 

huge tuning parameters Xh the estimator will return the estimate f3f^ = 0 = (0,..., 0)' where all parameters are 
shrinked to zero. In general the lasso estimator is a biased estimator. However, it can be tuned in such a way that 
the out-of-sample performance is better than the OLS estimator performance. The lasso estimator has an automatic 
sparsity property, so many of the potential parameters are not included in the finally estimated model. The larger A^ 
the less parameters are included in the estimated model. But lasso does not only sparsify, it also regularizes the model. 
So the estimated parameter that are not estimated to be exactly zero are usually smaller (in absolute values) than 
the corresponding OLS solution with the same parameters. To be precise: the in-sample residual sum of square (RSS) 
is larger for the lasso solution than for the corresponding OLS solution with the same non-zero parameters. From a 
Baysian point of view lasso acts as a prior on the parameters that is centered at zero. For more information on the 
lasso algorithm from the theoretical and applied point of view we recommend the recently published lasso overview 
book Hastie et al. (2015). 

As estimation algorithm the very fast coordinate descent algorithm of Friedman et al. (2007) is available. It solves 
the lasso problem on a given grid A of A-values that is commonly chosen as a exponential grid. For this exponential 
grid A it is important that the corresponding solutions cover the full range from the P — 0 solution to a solution that 
contains all possible parameters as the OLS solution. In detail we define the grid as A = {2^\k G G}U {0} where G is 
an equidistant grid from 4 to -15 of length 50. We use the R package glmnet for implementation, see e.g. Friedman 
et al. (2010) for more details. 

We choose the tuning parameters A/i G A for each hour h by minimizing the Bayesian information criterion (BIG) 
for each of the 24 single models. The BIG is a very conservative information criteria that avoids overestimation. 
However other information criteria, like the Akaike information criterion (AIG) or cross-validation based criteria, are 
possible as well. 

3 Benchmarks models 

As benchmarks we allow only those processes that model linear dependency structure in the data. These models 
are usually those ones that can be estimated and forecasted in a very fast way, like the proposed model with the 
corresponding estimation method. Furthermore, we take suitable representatives of the introduced two model classes. 
As we consider the lasso based model without (1) and with (2) day-of-the-week effects, we consider for the benchmark 
models a similar weekday modeling extension to keep comparability. 

3.1 24 univariate AR(p) models 

The model of 24 univariate AR(p) models is given by 

p 

Pd,h = 4>h,0 + ^ 4'h,kiPd-k,h — </>/i,o) + £d,h (7) 

k=l 

where {ed,h)d£Z is i.i.d. with Sd^h ~ -^(OjCr^)- It models the linear autoregressive dependence structure of each price 
to the corresponding hour well. We estimate the model by solving the Yule-Walker equations, which solution is 
guaranteed to be stationary. We select the model orders ph by minimizing the Akaike information criteria (AIG) given 
a maximal upper bound of ph by Pmax = 50. For the model version with weekday effects we consider the equation 

6 p / 6 \ 

Pd,h —^h,Q T ^ ^ Hj' (d) -t- ^ ^ kl^h,k I Pd—k,h ^ ^ J Hj (d ^) (8) 

i=0 k=l y j=0 j 

We estimate model (8) in a two-step approach. First we estimate the day-of-the-week coefficients iph,j by OLS and 
solve the Yule-Walker equations for the resulting residuals. As for model (7) we minimize the AIG with an upper 
bound for the order ph of Pmax = 50. Note that model (8) is just one possibility to include the day-of-the-week effects. 
Other ones are possible as well, but give similar results. We consider the chosen option as it allows for the application 
of the Yule-Walker equations as equivalently done for (7). 

3.2 Autoregressive expert models 

As mentioned in many applications, like in Weron and Misiorek (2008) or Maciejowska et al. (2015), simple autoregres¬ 
sive models with a fixed regression structure are used. As the assumed regression structure is usually build on some 
prior knowledge of experts, we refer these models as expert models. We consider two expert models as benchmarks, 
a simple expert model and a weekday based expert model. They are given by 


Pd,h = 4'h,0 + (t>h,lPd-l,h + 4’h,2Pd-2,h + 4>h,7Pd-7,h + Sd,h, 

Pd,h = 4'h,0 + 4‘h,lPd-l,h + (t>h,2Pd-2,h + <Ph,7Pd-7,h + '0oIIA)('^) + '4’lWi{d) -|- pQWe{d) + Sd,h- 


(9) 

( 10 ) 
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The latter one matches the first one but adds linear weekday effects on Sunday, Monday and Saturday. Here we 
assume for the error process again that {ed,h)d& is i-i-d. with Sd^h -^(0, cr^)- We estimate the parameters by solving 
the corresponding least square problem. The chosen model of Weron and Misiorek (2008) is just a representative of an 
expert model. There are several other versions available as well, they tend to slightly differ in some model components. 

For example Karakatsani and Bunn (2008a) propose a model with a linear price dependency on a daily lag of order 
1 and 7, but not at lag 2. Note that model (9) is nested in model (7), which is nested in the proposed model (1). 
Thus, we expect similar behavior, if the lags of the nested models capture the relevant features. These models with 
a simpler structure have the advantage of being easier to interpret and implement. Furthermore, the overestimation 
issue is negligible. 

3.3 Univariate AR(p) models 

This is a model of the second model class. Ziel et al. (2015a) and Ziel et al. (2015b) show that the univariate AR(p) 
model can be a powerful forecast method. It performs often even better than several modern models. Denote Pt = Pd,h 
with t = 24d + j for j € {0,..., 23} a univariate time series representation of electricity price at day d and hour h. 
The considered model and its weekday effect version are given by 

p 

Pt =(t>0 + 4>k{Pt-k — 4 * 0 ) + £t, 

k=l 

6 _ P / 6 _ \ 

Pt =(t >0 + X! + X! ^'kPt-k - X! V'itF, (J(t) - k)-(j)o\ +et 

j=0 k=l y j-O J 

where {ed,h)dei- is i-i-d. with ed,h ~ -^(0, cr^) and d{t) = LV24J gives the day d that corresponds to time point t. Again 
we estimate the model (11) by solving the Yule-Walker equations and minimize the AIC with respect to a maximal 
order of Pmax = 700. Note that a maximal oder larger than 168 is required to capture any weekly dependency structure. 

As for model (8) we estimate the weekday version (12) by the two-step approach. We solve the OLS regression for the 
day-of-the-week parameters and solve the Yule-Walker equations of the residuals by minimizing the AIC with maximal 
order of Pmax = 700. 

3.4 PCA based autoregressive models 

The model with 24 univariate AR(p) models is a special case of the 24-dimensional AR(p) process given by 

p 

-Pd = 00 + ^ ^k{Pd-k — 0o) + (13) 

with 24-dimensional intercept vector 0q, 24 x 24-dimensional autoregressive parameter matrices Pd = (Pd, 0 ; - - - j Pd, 23 )' 
and Ed N 2 i{ 0 , S). Such a general 24-dimensional AR(p)-process has 24-|- 24 x 24 x p = 24-|- 576p parameters. This 
is usually quite a lot and leads quickly to models with too many parameters. Note that model (13) is able to capture 
the lagged cross-hour dependency structure as seen in Figure 1. The proposed model (1) is nested in (13) (and vice 
versa) with an appropriate choice of I{h, 1). As the lasso estimation technique overcomes the overestimation problem 
the same features can be covered. Another approach that deals with this over-parametrization problem is principal 
component analysis (PCA) based. 

The approach is based on a dimension reduction of the 24-dimensional process by PCA. Maciejowska and Weron 
(2013) and Maciejowska and Weron (2015) use it and a generalization to forecast the UK APX electricity spot prices. 
They show that the PCA approach can be suitable to capture the intraday relationship between the electricity prices. 
The principal component decomposition of Pd is given by 

Pd = PFd (14) 

with r as 24X 24-dimensional loading matrix and Fd as factor scores at day d. As usual in PCA analysis we assume 
that without loss of generality the columns of P are ordered by the absolute value of their corresponding eigenvalue. 

For the dimension reduction we consider only the first K principal components of the scores which we denote by FK,d- 
We denote the corresponding AT x 24-dimensional reduced loading matrix by Fjf. 

For forecasting the scores we employ a multivariate AR (VAR) model to the AT-dimensional vector FK,d- Thus, 
we assume that the factor scores follow a process that is given by 

p 

FK,d = (t>K ,0 + ^ ^K,k{FK,d-k — (t>K,o) + ^K,d (15) 

k=l 

with AT-dimensional intercept vector K x K dimensional autoregressive parameter matrices ^K,k and £K,d as 

AT-dimensional zero-mean error term vector. Now the dimension reduced model has only K + K x K x p parameters. 


( 11 ) 

( 12 ) 
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As it holds Pd = + J2k=i ^K^K, k{FK,d—k — 4^k,o) a ^K^K,d model (15) is a special case of model (13), but 

also of the proposed model (1) if I(/i, 1) is chosen correctly. 

Again we estimate model (15) by solving the multivariate Yule-Walker equations and minimizing the AIC with 
respect to a maximal order of Pmax = 50. We perform the estimation for K G {2,..., 12}. Given a forecast FK,d of the 
factor scores for day d and an estimate F of the loading matrix we can easily get the price forecast by Pd = TkFk^ 
using the plug-in principle. 

The day-of-the-week version of model (15) is given by 

6 P / 6 \ 

FK,d =(t>K.O + I ^K,kF K,d-k — — k) j +eK,d (16) 

j=o k=i y i=o J 

where are the AT-dimensional vectors of weekday parameter coefficient, and 0^ g is again the A'-dimensional 
intercept vector, ^K,k the K x K autoregressive parameter matrices and £K,d as AT-dimensional error process. As for 
models (8) and (12) we estimate the day-of-the-week model (12) by the two-step approach. First we solve the OLS 
problem for the weekday parameters and take the Yule-Walker equations solution of the residuals by minimizing the 
AIC with maximal possible order of Pmax = 50 and for AT € {2,..., 12}. 


4 Empirical study and results 

We perform an out-of-sample forecast study by using a rolling window. As in-sample period we select always D = 
730 = 2 X 365 days (usually two years) of data for the estimation. Given the estimated model we perform an out-of- 
sample forecast for the next 24 values. The considered data range is from 17.12.2009 to 12.08.2014. Thus, the full 
out-of-sample period is from 17.12.2011 to 12.08.2014 which is exactly N = 969 days. 

As forecasting measure we consider the overall MAE and RMSE, as they are also suitable on markets with negative 
prices. The MAE is a more robust measure, whereas the RMSE is the optimal measure for least square problems and 
more sensitive to outliers. They are given by 


MAE = 


1 

2^ 


N 23 

\PD+n,h - 


n—1h—Q 


P, 


Z)+n,/i 


RMSE = 




1 

24/V 


N 23 

Y Y lPD+n,h 

n—1h—0 


PD+n,hl^ 


(17) 

(18) 


where PD+n,h denotes the forecasted price of day D + n. However, to better understand how the models perform over 
the day, we use the hourly MAE and RMSE at hour h, denoted by MAE?i and RMSE^^. They are defined through 


MAE^ = 


N 


N 


Y\^^ 


D-\-n,h 


-P 


D-\-n,h \ 5 


RMSE,, 




1 

N 


N 

Y, \PD+n,h — PD+n,h\'^- 

n—1 


(19) 

( 20 ) 


We conduct the out-of-sample study for a selected number of electricity markets and all considered models. For the 
PGA based models in (15) we only report the model with minimal MAE for each market. The selected electricity 
prices are the German/Austrian EXAA price, the German/Austrian and Swiss EPEX prices, the Belgian BELPEX 
price, the Dutch APX price, the Danish West, Danish East and Sweden(4) Nordpool prices, the Polish POLPX price 
and the Czech OTE price. A summary of the considered prices with the subsequently used abbreviations is given in 
Table 1. Note that all prices are hourly day-ahead electricity prices. For comparison issues we converted all prices 
in EUR/MWh. If the reported currency is not EUR we consider the corresponding end of day exchange rate of the 
local currency to convert the prices to EUR/MWh. A selection of different electricity prices is suitable to check if 
the intraday structure on different electricity markets is similarly or not. It also acts as a robustness check for the 
proposed lasso based estimation method. 

For the PCA-based models the optimal number of factors K are between 3 and 10. They are also listed in Table 
1. Note that Maciejowska and Weron (2015) get an optimal PGA-based model with 3 principal components for the 
half-hourly volume-weighted APX UK spot log-price. They also provide a nice interpretation for the factors. 

Subsequently, we denote model (1) by lasso, (7) by 24d.AR, (9) by exp.AR, (11) by AR(p) and (15) by PGA*. 
For the day-of-the-week model version we add the suffix ’-wd’ to the model abbreviation. The computed MAE and 
RMSE values are given in Table 2. There we also report the estimated standard deviations of MAE and RMSE which 
are computed using residual based bootstrap with a bootstrap sample size of 10000. These standard deviations can 
be used as significance criterion as in Ziel et al. (2015b). For example if the MAE and RMSE would be normally 
distributed, then the 1.96-sigma range around the mean matches the symmetric 95% confidence interval. Even though 
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Exchange 

Region 

Abbreviation 

PCA* 

PCA*-' 

EXAA 

Germany&Austria 

EXAA.DE&AT 

7 

8 

EPEX 

Germany&Austria 

EPEX.DE&AT 

6 

10 

EPEX 

Switzerland 

EPEX.GH 

8 

6 

BELPEX 

Belgium 

belpex.be 

6 

6 

APX 

Netherlands 

APX.NL 

6 

8 

Nordpool 

Denmark (West) 

NP.DK.West 

7 

3 

Nordpool 

Denmark (East) 

NP.DK.East 

4 

4 

Nordpool 

Sweden (4) 

NP.SW 

6 

4 

POLPX 

Poland (Auction I) 

polpx.pl 

8 

8 

OTE 

Gzech Republic 

OTE.GZ 

7 

9 


Table 1: Summary of the considered electricity markets with abbreviations and optimal number of factors for the PCA 
based models. 


the distribution is likely not normal we can use the Chebyshev inequality or extensions to make non-parametric 
statements about the confidence interval. In Table 2 we underlined all models that are not significantly worse than 
the best model, where we take the 2-sigma range of the best model as significance indicator. 

In Table 2 we observe that concerning the MAE the lasso based approaches perform very well. The proposed model 
(lasso-wd) in equation (2) with weekdays effects is always the best model, except for the BELPEX and POLPX case. 
But even there it not significantly worse than the best models which are the lasso without weekday effects for BELPEX 
and the 24-dimensional AR model with weekday effects for the POLPX. It is remarkable that model (2) with weekdays 
effects is sometimes significantly better than all other considered models. In several cases (EPEX.AT&AT, EPEX.CH, 
APX.NL, NP.DK.East, OTE.CZ) there is only one model that is not significantly worse. In contrast, concerning the 
RMSE values the picture is not that clear. Still, the lasso-wd model (2) is most of the time the best model. However, 
as mentioned the RMSE is sensitive to outliers. So for example the Nordpool Denmark West price has 5 extreme 
prices between 1900EUR/MWh and 2000EUR/MWh on June 7, 2013. This leads to a drastic increase in the RMSE, 
as no (considered) model is able to capture these extreme price spikes. Concerning the RMSE the AR(p)-wd seems 
to perform relatively well and scores best for 3 of the 10 considered markets. The reason might be that the AR(p)-wd 
has a relatively low outlier sensitivity, as it is a member of the second model class which contains big single equation 
models and use all observations as input. Having a lot of observations helps these models to learn that the outliers 
are a rare event. In contrast, all models of the first class that use 24 small models (as the proposed lasso) have the 
problem that an outlier gets much more attention as less observations are available. 

For the MAE?, and RMSE?, we discuss detailed results only for the German/Austrian EPEX spot price, which is 
the largest day-ahead electricity market in Europe (regarding the traded volume). The MAE?, and RMSE?, results are 
given in Figure 3. Results for the other markets are similar. The corresponding graphs are given in the supplementary 
material. There we also present the correlation plots and the weekly mean graphs as in Figures 1 and 2 for all markets. 


-e- 

lasso . - 

exp.AR —0— 

PCA* 

-A- 

lasso-wd 

exp.AR-wd — e— 

PCA*-wd 


24d.AR -B- 

AR(p) 



24d.AR—wd — 

AR(p)-wd 




0 2 4 6 8 10 13 16 19 22 

(a) MAE,, 



0 2 4 6 8 10 13 16 19 22 


(b) RMSE,, 


Figure 3: Hourly performance measures for EPEX spot price of Germany and Austria. 
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MAE 

lasso 

lasso-wd 

24d.AR 

24d.AR-wd 

exp.AR 

exp.AR-wd 

AR(p) 

AR(p)-wd 

PCA* 

PCA*-wd 

EXAA.DE&AT 

4.89 

4.61 

5.30 

4.96 

6.09 

5.27 

4.87 

4.68 

5.11 

4.72 

(.033) 

(.031) 

(.035) 

(.034) 

(.041) 

(.035) 

(.033) 

(.032) 

(.034) 

(.032) 

EPEX.DE&AT 

6.46 

6.01 

6.73 

6.32 

7.43 

6.69 

6.39 

6.17 

6.62 

6.09 

(.050) 

(.048) 

(.054) 

(.052) 

(.058) 

(.054) 

(.050) 

(.049) 

(.052) 

(.050) 

EPEX.CH 

5.44 

5.15 

5.82 

5.36 

6.68 

5.65 

5.37 

5.23 

5.93 

5.38 

(.047) 

(.045) 

(.048) 

(.045) 

(.054) 

(.047) 

(.045) 

(.043) 

(.052) 

(.048) 

belpex.be 

6.75 

6.78 

6.91 

7.19 

7.26 

7.38 

7.21 

7.48 

7.05 

7.27 

(.052) 

(.052) 

(.053) 

(.054) 

(.056) 

(.056) 

(.053) 

(.053) 

(.052) 

(.054) 

apx.nl 

5.11 

4.95 

5.26 

5.04 

5.62 

5.30 

5.13 

5.01 

5.25 

5.03 

(.034) 

(.033) 

(.035) 

(.033) 

(.038) 

(.035) 

(.034) 

(.033) 

(.034) 

(.033) 

NP.DK.West 

6.89 

6.62 

7.21 

6.95 

8.03 

7.86 

7.61 

7.42 

6.83 

6.78 

(.193) 

(.191) 

(.187) 

(.186) 

(.298) 

(.328) 

(.194) 

(.194) 

(.188) 

(.187) 

NP.DK.East 

6.31 

6.01 

6.42 

6.18 

6.66 

6.40 

6.13 

6.07 

6.50 

6.35 

(.069) 

(.067) 

(.055) 

(.054) 

(.057) 

(.057) 

(.053) 

(.052) 

(.056) 

(.055) 

NP.SW4 

3.28 

3.20 

3.67 

3.60 

3.73 

3.54 

3.48 

3.52 

3.78 

3.74 

(.038) 

(.038) 

(.039) 

(.037) 

(.039) 

(.038) 

(.039) 

(.038) 

(.039) 

(.039) 

polpx.pl 

2.78 

2.68 

2.82 

2.65 

3.19 

2.73 

2.83 

2.77 

3.04 

2.81 

(.030) 

(.029) 

(.027) 

(.027) 

(.030) 

(.028) 

(.027) 

(.027) 

(.030) 

(.028) 

OTE.cz 

5.49 

5.10 

5.67 

5.26 

6.52 

5.63 

5.34 

5.11 

5.60 

5.20 

(.039) 

(.036) 

(.041) 

(.039) 

(.047) 

(.042) 

(.038) 

(.037) 

(.040) 

(.038) 

RMSE 

lasso 

lasso-wd 

24d.AR 

24d.AR-wd 

exp.AR 

exp.AR-wd 

AR(p) 

AR(p)-wd 

PCA* 

PCA*-wd 

EXAA.DE&AT 

7.02 

6.63 

7.52 

7.12 

8.70 

7.58 

7.01 

6.73 

7.33 

6.79 

(.070) 

(.067) 

(.068) 

(.066) 

(.081) 

(.068) 

(.067) 

(.066) 

(.072) 

(.070) 

EPEX.DE&AT 

10.06 

9.47 

10.70 

10.24 

11.63 

10.71 

10.02 

9.70 

10.45 

9.82 

(.275) 

(.282) 

(.308) 

(.323) 

(.284) 

(.301) 

(.265) 

(.272) 

(.296) 

(.324) 

EPEX.CH 

9.03 

8.60 

9.35 

8.72 

10.58 

9.16 

8.77 

8.48 

9.93 

9.18 

(.220) 

(.227) 

(.211) 

(.211) 

(.212) 

(.201) 

(.208) 

(.210) 

(.260) 

(.267) 

belpex.be 

10.33 

10.37 

10.58 

10.98 

11.21 

11.30 

10.79 

11.02 

10.58 

10.90 

(.239) 

(.239) 

(.234) 

(.243) 

(.232) 

(.231) 

(.229) 

(.222) 

(.210) 

(.215) 

apx.nl 

7.26 

7.02 

7.44 

7.14 

8.00 

7.50 

7.26 

7.09 

7.40 

7.09 

(.079) 

(.074) 

(.077) 

(.073) 

(.087) 

(.076) 

(.078) 

(.076) 

(.081) 

(.079) 

NP.DK.West 

30.21 

29.94 

29.47 

29.37 

46.27 

50.67 

30.69 

30.50 

29.61 

29.35 

(5.845) 

(5.884) 

(5.934) 

(5.956) 

(7.033) 

(8.386) 

(5.662) 

(5.697) 

(6.004) 

(6.016) 

NP.DK.East 

12.38 

11.92 

10.59 

10.26 

11.07 

10.83 

10.34 

10.06 

10.77 

10.52 

(.917) 

(.826) 

(.285) 

(.297) 

(.281) 

(.289) 

(.311) 

(.311) 

(.288) 

(.297) 

NP.SW4 

6.76 

6.64 

6.97 

6.72 

7.06 

6.82 

6.90 

6.78 

7.08 

7.10 

(.223) 

(.226) 

(.225) 

(.229) 

(.217) 

(.224) 

(.241) 

(.243) 

(.220) 

(.221) 

polpx.pl 

5.37 

5.17 

5.06 

4.94 

5.61 

5.08 

4.98 

4.92 

5.52 

5.19 

(.156) 

(.145) 

(.127) 

(.130) 

(.133) 

(.133) 

(.114) 

(.115) 

(.139) 

(.138) 

OTE.cz 

8.08 

7.57 

8.46 

7.98 

9.66 

8.54 

7.90 

7.58 

8.25 

7.81 

(.140) 

(.137) 

(.155) 

(.154) 

(.157) 

(.145) 

(.121) 

(.123) 

(.136) 

(.150) 


Table 2: MAE and RMSE values for the considered markets and models. The corresponding estimated standard 
deviations are given in parenthesis. The best model is highlighted in bold font. All models that are not significantly 
worse than best (indicated by the 2-sigma range of the best model) are underlined. 


In Eigure 3 we observe that the daily pattern for the MAE/j and RMSE^ are similar. In general we see that 
the model with weekday effects estimated by lasso performs as one of the best models for all hours over the day. It 
seems that during the morning hours it outperforms the other models most clearly. It is interesting that in the night 
hours from 0:00 to 5:00 the lasso based models perform similarly to the univariate AR(p). But this matches well the 
observation from the beginning that these hours are highly correlated with the last hours of the previous day. The 
univariate AR(p) and the lasso models can capture this behavior. Moreover, in the evening hours (esp. 21:00 to 23:00) 
the 24-dimensional AR model and the 24-dimensional autoregressive expert model perform as good as the lasso based 
models. Again, this coincides with the observation from Figure 1, that these hours are strongly correlated with the 
same hours of the previous days. The mentioned models are all able to cover this dependency well. Furthermore, we 
see a relatively well performance of the PCA*-wd model. We know from Table 2 that for the German/Austrian EPEX 
price only the PCA*-wd model performs not significantly worse. The behavior is similar as both core models (ignoring 
the weekday effects) are nested in (13). However, regarding the cross-country Table 2, we know that in most cases 
the lasso-wd model outperforms the PCA*-wd model. The reason is likely that cross-hour dependency is modeled 
directly instead of using the reduced PCA representation. It might happen that relevant information gets lost due to 
the dimension reduction technique. 

Another interesting fact in Figure 3 is that it seems that in the morning hours around 6:00 to 7:00 the lasso-wd 
performs much better in the MAE^, score than in the RMSE/,. However, having a closer look at the graph it turns out 
the AR(p) process catches up in terms of RMSE/,. The reason might be again that the AR(p) based models are less 
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outlier sensitive than the other considered models. Interestingly this pattern does not hold for all considered markets, 
but we observe the same fact for EXAA.DE+AT, EPEX.CH, APX.NL, NP.DK.East and OTE.CZ. 

Furthermore, we can see in Figure 3 that it seems that the weekday dummies in the model (2) improve the 
performance significantly compared to model (1). We will understand this better when we analyze the intraday 
structure and the variable importance. 

Hence, additionally to the out-of-sample study we have a closer look at the intraday structure of the electricity 
prices by using the proposed model (2). As the lasso estimator sparsities and regularizes, we can study the finally 
selected parameters and their magnitude. For a measure of importance of a parameter f3h,i resp. /3h,i we consider ih,i, 
the fraction of the absolute standardized parameter of a model to the sum of all absolute standardized parameters: 






( 21 ) 


Obviously, we have that 0 < Lh,i < 1 and ‘'h,i = 1- The larger Lh.,i the larger the relative impact to the electricity 
price of hour h. We estimate ih,i for model (2) which includes the weekday dummies for all considered markets by the 

corresponding plug-in estimator using the estimated standardized parameter vector /3^. Given the estimated 
values, we can sort them decreasingly. Then the largest Th,i value matches the most relevant parameter in the model 
for hour h. 

Again, we only discuss the results for the German/Austrian EPEX spot price in detail. In Table 3 we present 
the 7 most important variables for each model of hour h and their Th,i value. The corresponding tables for the other 
electricity prices are given in the supplementary material. In Table 3 we can observe an interesting pattern that we 



Importance: I 

Importance: 2 

Importance: 3 

Importance: 4 

Importance: 5 

Importance: 6 

Importance: 7 

0 

23@I (70.6) 

19@I ( 5.9) 

I@I ( 5.7) 

18@I ( 4.1) 

21@1 ( 3.6) 

4@1 ( 3.1) 

Sat 

(2.2) 

1 

23@I (67.3) 

19@I ( 8.6) 

3@I ( 7.1) 

4@1 ( 5.9) 

Sun ( 4.7) 

5@1 ( 2.6) 

6@1 

(2.0) 

2 

23@I (56.3) 

4@I ( 9.4) 

19@I ( 8.9) 

3@1 ( 7.9) 

Sun ( 6.6) 

22@7 ( 5.1) 

2@25 

( 1-6) 

3 

23@I (45.0) 

4@I (14.3) 

3@I (12.2) 

19@1 ( 8.5) 

Sun ( 7.3) 

22@7 ( 4.3) 

3@6 

(2.5) 

4 

23@I (44.8) 

4@I (13.8) 

19@I (10.3) 

Sun (10.2) 

3@1 ( 7.4) 

22@7 ( 2.5) 

4@25 

(2.4) 

5 

23@I (34.1) 

Sun (15.5) 

5@I ( 9.7) 

4@1 ( 8.2) 

Sat ( 7.4) 

19@1 ( 7.2) 

22@7 

(3.8) 

6 

Sun (20.4) 

23@I (20.3) 

Sat (12.8) 

6@1 ( 9.4) 

Mon ( 8.7) 

20@1 ( 6.0) 

23@2 

(4.1) 

7 

Sun (23.9) 

Sat (15.9) 

23@I (I2.I) 

Mon ( 7.8) 

20@1 ( 6.3) 

6@1 ( 5.3) 

I9@l 

(4.5) 

8 

Sun (28.3) 

Sat (15.9) 

20@I ( 8.5) 

23@1 ( 7.9) 

18@1 ( 7.0) 

Mon ( 5.4) 

21@1 

(5.2) 

9 

Sun (25.9) 

20@I (13.3) 

Sat (13.3) 

23@1 ( 7.0) 

22@1 ( 6.3) 

Mon ( 6.1) 

I8@l 

(5.5) 

10 

Sun (21.8) 

Sat (11.9) 

22@I ( 9.6) 

20@1 ( 9.0) 

Mon ( 6.6) 

13@1 ( 6.1) 

23@1 

(5.6) 

11 

Sun (17.2) 

22@1 (12.3) 

Sat (11.7) 

13@1 ( 9.2) 

Mon ( 8.0) 

23@1 ( 5.8) 

Il@7 

(4.1) 

12 

Sun (15.6) 

13@1 (12.0) 

Sat (10.7) 

22@1 (10.2) 

Mon ( 8.1) 

23@1 ( 8.0) 

I7@l 

(4.9) 

13 

Sun (18.6) 

Sat (11.2) 

13@1 (11.1) 

22@1 (10.0) 

Mon ( 7.9) 

23@1 ( 7.3) 

I7@l 

(6.8) 

14 

Sun (18.7) 

Sat (10.1) 

Mon ( 8.2) 

23@1 ( 8.2) 

16@1 ( 8.1) 

22@1 ( 7.6) 

I3@l 

(4.2) 

15 

Sun (19.8) 

16@1 (12.5) 

Sat ( 8.9) 

Mon ( 8.0) 

23@1 ( 7.1) 

17@6 ( 6.0) 

I7@l 

(4.9) 

16 

Sun (19.3) 

17@1 (12.5) 

Sat ( 8.8) 

16@1 ( 8.5) 

Mon ( 8.4) 

23@1 ( 7.1) 

I7@6 

(6.7) 

17 

17@1 (26.4) 

Sun (16.1) 

Mon (10.1) 

17@6 (10.0) 

Sat ( 6.7) 

17@7 ( 6.3) 

I8@l 

(6.0) 

18 

18@1 (33.8) 

Sun (14.1) 

Mon ( 9.9) 

18@6 ( 9.4) 

Sat ( 6.0) 

18@2 ( 5.9) 

I8@5 

(4.6) 

19 

19@1 (34.0) 

Mon ( 9.4) 

19@6 ( 8.5) 

Sun ( 8.2) 

19@2 ( 7.6) 

8@6 ( 5.0) 

I9@7 

(4.7) 

20 

20@1 (23.9) 

20@2 (10.4) 

Sat ( 9.5) 

Mon ( 7.8) 

21@1 ( 7.5) 

20@6 ( 7.5) 

Sun 

(7.0) 

21 

21@1 (22.8) 

Sat (10.9) 

21@2 ( 9.9) 

Mon ( 6.2) 

21@28 ( 4.7) 

Sun ( 4.6) 

23@1 

(4.6) 

22 

22@1 (30.6) 

22@2 ( 9.5) 

Sat ( 7.7) 

23@3 ( 6.5) 

22@28 ( 6.1) 

22@5 ( 5.1) 

22@35 

(5.0) 

23 

23@1 (24.8) 

22@1 (17.5) 

23@3 (10.7) 

22@2 ( 9.6) 

Sat ( 5.7) 

0@5 ( 5.0) 

I2@7 

(4.4) 


Table 3: Most relevant coefficients in model (2) for each hour of the German/Austrian EPEX price. H@L represents 
the estimated parameter for Yd-L,H for a model on Yd^h- Sun, Mon, ..., Sat represent the weekday dummies. The 
corresponding Thd value is given in parenthesis. 


have similarly seen already in the introduction. The most relevant parameters for the first 6 hours of the day is the 
last observed value of the previous day (which corresponds to Pd- 1 , 23 )- For the first hour this importance is very 
distinct. The last observation Pd- 1,23 of the previous day contributes 70.6% to the full model. From the physical point 
of view this makes sense. Obviously, the physical delivery of the auction result for a specific hour is based on related 
processes like the electricity load or wind energy production. These situations should not differ that much between 
23:00 and the next day 0:00. So a large impact of the 23:00 price to the price of the next day at hour 0:00 is plausible. 
This observation seems to be very stable across all countries. So we observe for all considered markets that the most 
relevant variable for the first 0:00 hour is always the last hour of the previous day. For all markets, except for the 
BELPEX.be, this is even the case for the first 3 hours (0:00, 1:00 and 2:00) of the day. A similar result is recently 
established in load forecasting by Dudek (2016). However, for electricity price forecasting this is not that obvious as 
we observe all 24 values at once. The last observation of the previous day is not more recent than the other ones. 
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The fact that for the first hours the last values of the previous day are very important underlines that the univariate 
AR(p) and the lasso based models perform well in this time horizon. 

For the last 7 hours of the day (from 17:00 to 23:00) the most relevant lag is always the value of the previous 
day at the same hour (which corresponds to Pd-i,h)- This matches well the observed fact that the 24-dimensional 
autoregressive models have a good out-of-sample performance for these hours. 

Moreover, we see that for the morning and afternoon hours from 6:00 to 16:00 the most important parameter is 
the weekday dummy on Sunday. So the weekday effect seems to be more important than the autoregressive impact. 
In general we observe that the weekday dummies have a larger impact during the working hours than in the night and 
evening hours. Thus, we can deduce that the weekly seasonal structure of the electricity price is more distinct in the 
working hours. This coincides with Figure 2a, where we see that the difference in the weekly mean is smaller during 
night and evening than the rest of the day. 

Additionally, we observe that only the weekday dummies for Sunday, Monday and Saturday occur in Table 3. Here 
it is interesting that the considered expert model from Weron and Misiorek (2008) uses exactly the same dummies. 
However, regarding the weekly sample means in Figure 2a this is not obvious. From this graph we would expect that 
the impact of Monday dummies is similar to Friday dummies. Furthermore the Monday dummies appear only for the 
hours from 6:00 to 21:00 which is definitively not expectable from Figure 2a which suggests a Monday dummy only 
from 0:00 to 5:00. 

Another fact that we can derive from Table 3 concerns the weekly lagged impact. In many models, e.g. the expert 
model from Weron and Misiorek (2008) a linear dependence on the weekly lagged electricity price is considered Pd- 7 ,h- 
However, in Table 3 it never appears in the top 4 for any hour. Only for some hours it is estimated to be the 5th, 
6th or 7th most important variable. This is an indication that given a proper short-term autoregressive and weekday 
modeling the modeled electricity price dependence from the weekly lagged price gets weaker. The last fact from Table 
3 that we want to mention concerns the evening hours. We observe that only for the evening hours from 18:00 to 23:00 
a dependence of electricity price of the day before yesterday (lag = 2) appears. Again, this dependency is always to 
the price at the same hour. 

5 Summary and conclusion 

We present a day-ahead electricity spot price model approach for hourly data that is able to capture the intraday 
dependency structure well, especially the time-varying cross-hour dependencies. Even though the considered regression 
model consists of 24 simple linear regression models it explains the variety in the price data. This relies heavily on the 
utilized lasso estimation technique. It automatically shrinks and sparsifies the model parameters. Thus, the estimated 
24 single models can vary over the day and capture the dependency in the data. 

We show that the dependency structure significantly changes over the days. For example in the German/Austrian 
EPEX market, the night hours from 0:00 to 5:00 have a strong dependency to the last price at 23:00 of the previous 
day. This result holds similarly across all considered European markets. In contrast, the price at the evening hours 
of from 18:00 to 23:00 of the German/Austrian EPEX market depends more on the price of the previous day at the 
same hour. We highlight that the majority of available electricity price model approaches are not designed to cover 
both features, whereas the lasso estimation technique captures both features in an efficient way. 

The conducted out-of-sample forecasting study shows that the considered lasso based approach performs strong in 
terms of forecasting accuracy. This is remarkable, if we have in mind that the estimation and forecasting procedure is 
very fast. Hence, the model is a suitable candidate for forecast combination approaches as in Bordignon et al. (2013) 
or Maciejowska et al. (2015). As the out-of-sample study shows similar performance results for all 10 considered 
electricity markets, it underlines that the flexible model approach can cover the local market specific dependency 
behavior. 

Future research can go in different directions. The likely most important issue is the covering of the non-linear 
behavior, especially those effects that induce price spikes. This should be combined with an approach that takes other 
sources like renewable energy feed-in or load into account. Especially the impact of wind power is very important 
to cover price spikes, see e.g. Keles et al. (2012). Furthermore, there should be more research done for probabilistic 
forecasting. 


References 

Abedinia, O., Amjady, N., Shafie-khah, M., and Gatalao, J. (2015). Electricity price forecast using combinatorial 
neural network trained by a new stochastic search method. Energy Conversion and Management^ 105:642-654. 

Bordignon, S., Bunn, D. W., Lisi, F., and Nan, F. (2013). Gombining day-ahead forecasts for British electricity prices. 
Energy Economics, 35:88-103. 

Gatalao, J., Pousinho, H., and Mendes, V. (2011). Hybrid wavelet-pso-anfis approach for short-term electricity prices 
forecasting. IEEE Transactions on Power Systems, 26(1):137-144. 



References 


12 


Chaabane, N. (2014). A hybrid arfima and neural network model for electricity price prediction. International Journal 
of Electrical Power & Energy Systems, 55:187-194. 

Conejo, A. J., Plazas, M., Espinola, R., Molina, A. B., et al. (2005). Day-ahead electricity price forecasting using the 
wavelet transform and ARIMA models. Power Systems, IEEE Transactions on, 20(2):1035-1042. 

Cruz, A., Munoz, A., Zamora, J. L., and Espinola, R. (2011). The effect of wind generation and weekday on Spanish 
electricity spot price forecasting. Electric Power Systems Research, 81(10):1924-1935. 

Dudek, G. (2016). Pattern-based local linear regression models for short-term load forecasting. Electric Power Systems 
Research, 130:139-147. 

Friedman, J., Hastie, T., Hofling, H., and Tibshirani, R. (2007). Pathwise coordinate optimization. The Annals of 
Applied Statistics, l(2):302-332. 

Friedman, J., Hastie, T., and Tibshirani, R. (2010). Regularization paths for generalized linear models via coordinate 
descent. Journal of statistical software, 33(1):1. 

Haldrup, N., Nielsen, F. S., and Nielsen, M. O. (2010). A vector autoregressive model for electricity prices subject to 
long memory and regime switching. Energy Economics, 32(5):1044-1058. 

Hastie, T., Tibshirani, R., and Wainwright, M. (2015). Statistical Learning with Sparsity: The Lasso and Generaliza¬ 
tions. CRC Press. 

Hu, L. and Taylor, G. (2014). A novel hybrid technique for short-term electricity price forecasting in uk electricity 
markets. Journal of International Council on Electrical Engineering, 4(2):114-120. 

Karakatsani, N. V. and Bunn, D. W. (2008a). Forecasting electricity prices: The impact of fundamentals and time- 
varying coefficients. International Journal of Forecasting, 24(4):764-785. 

Karakatsani, N. V. and Bunn, D. W. (2008b). Intra-day and regime-switching dynamics in electricity price formation. 
Energy Economics, 30(4):1776-1797. 

Keles, D., Genoese, M., Most, D., and Fichtner, W. (2012). Gomparison of extended mean-reversion and time series 
models for electricity spot price simulation considering negative prices. Energy Economics, 34(4):1012-1032. 

Ludwig, N., Feuerriegel, S., and Neumann, D. (2015). Putting big data analytics to work: Feature selection for 
forecasting electricity prices using the lasso and random forests. Journal of Decision Systems, 24(l):19-36. 

Maciejowska, K., Nowotarski, J., and Weron, R. (2015). Probabilistic forecasting of electricity spot prices using Factor 
Quantile Regression Averaging. International Journal of Forecasting. 

Maciejowska, K. and Weron, R. (2013). Forecasting of daily electricity spot prices by incorporating intra-day re¬ 
lationships: Evidence form the UK power market. In European Energy Market (EEM), 2013 10th International 
Conference on the, pages 1-5. IEEE. 

Maciejowska, K. and Weron, R. (2015). Short-and mid-term forecasting of baseload electricity prices in the UK: 
The impact of intra-day price relationships and market fundamentals. Power Systems, IEEE Transactions on, 
PP(99):1-12. 

Nowotarski, J. and Weron, R. (2014). Gomputing electricity spot price prediction intervals using quantile regression 
and forecast averaging. Computational Statistics, pages 1-13. 

Paraschiv, F., Erni, D., and Pietsch, R. (2014). The impact of renewable energies on EEX day-ahead electricity prices. 
Energy Policy, 73:196-210. 

Paraschiv, F., Fleten, S.-E., and Schiirle, M. (2015). A spot-forward model for electricity prices with regime shifts. 
Energy Economics, 47:142-153. 

Shafie-Khah, M., Moghaddam, M. P., and Sheikh-El-Eslami, M. (2011). Price forecasting of day-ahead electricity 
markets using a hybrid forecast method. Energy Conversion and Management, 52(5):2165-2169. 

Singhal, D. and Swarup, K. (2011). Electricity price forecasting using artificial neural networks. International Journal 
of Electrical Power & Energy Systems, 33(3):550-555. 

Tibshirani, R. (1996). Regression shrinkage and selection via the lasso. Journal of the Royal Statistical Society: Series 
B (Statistical Methodology), 58(l):267-288. 

Weron, R. (2014). Electricity price forecasting: A review of the state-of-the-art with a look into the future. International 
Journal of Forecasting, 30(4):1030-1081. 



References 


13 


Weron, R. and Misiorek, A. (2008). Forecasting spot electricity prices: A comparison of parametric and semiparametric 
time series models. International Journal of Forecasting, 24(4):744-763. 

Ziel, F. (2015). Iteratively reweighted adaptive lasso for conditional heteroscedastic time series with applications to 
ar-arch type processes. Computational Statistics and Data Analysis (in press). 

Ziel, F., Steinert, R., and Husmann, S. (2015a). Efficient modeling and forecasting of electricity spot prices. Energy 
Economics, 47:98-111. 

Ziel, F., Steinert, R., and Husmann, S. (2015b). Forecasting day ahead electricity spot prices: The impact of the exaa 
to other european electricity markets. Energy Economics, 51:430-444. 



6 Figures 


14 


Appendix 

6 Figures 

6.1 Correlation plots 
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Figure 4: Sample correlations of Cor {Pd, h, Pd-i,i) 
12.08.2014. 


nnnnnnnnnnnnnnnn n nnnnn n! i 
□□□□□□□□□□□□□□□□□□□□□nnn 
□nnnnnnnnnnnnnnnnnnnnnnn 
nnnnnnnnnnnnnnnnnnnnnnnn 
nnnnnnnnnnnnnnnnnnnnnnnn 
nnnnnnnnnnnnnnnnnnnnnnnn 
nnnnnnnnnnnnnnnn n nnnnnnn 

□□□□□□□□□□□□□□□□□□□□nnnn 

nnnnnnnnnnnnnnnn n nnnnnnn 

□□□□□□□□□□□□□□□□□□□□nnnn 

□□□□□□□□□□nnnnnnnnnnnnnn 

nnnnnnnnnnnnnnnn n nnnnnnn 

□nnnnnnnnnnnnnnn n nnnnnnn 

□□□□□□□□□□nnnnnnnnnnnnnn 

nnnnnnnnnnnnnnnn n nnnnnnn 

□□□□□nnnnnnnnnnn n nnnnnnn 

□□□□□□□□□□nnnnnnnnnnnnnn 

nnnnnnnnnnnnnnnn n nnnnnnn 

□□□□□□□□□□nnnnnn n nnnnnnn 

□□□□□□□□□□nnnnnnnnnnnnnn 

nnnnnnnnnnnnnnnn n nnnnnnn 

□□□□□□□□□□□□□□□□□□□□nnnn 

□□□□□□□□□□□□□□□□□□□□nnnn 

□□□□□□□□□□□□□□□□□□□□nnnn 


(b) EPEX - Germany and Austria 


Pd- 


d-1,1 


1 

0.8 

0.6 

0.4 

0.2 

0 

- 0.2 

-0.4 

- 0.6 

- 0.8 

-r 




^(Ntxlxtln'Ot''00o^o—'(Nmxt- 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

•=11 
cvl2 
^ 13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 


□□□□□□□□□□□□□□□□□□□□□nnn 

□□□□□□□□□□□□□□□□□□□□□nnn 

□□□□□□□□□□□□□□□□□□□□□□□□ 

□□□□□□□□□□□□□□□□□□□□□□□□ 

□□□□□□□□□□□□□□□□□□□□□□□□ 

□□□□□□□□□□□□□□□□□□□□□□□□ 

nnnnnnnnnnnnnnnn n nnnnnnn 


□□□□nnnn 

□□□□nnnn 


□□□□□ 

□□□□□ 


□nnnnnnnnnnnnnnnnnnnnnnn 

□□□□□nnnnnnnnnnnnnnnnnnn 

□nnnnnnnnnnnnnnnnnnnnnnn 

□nnnnnnnnnnnnnnnnnnnnnnn 

□nnnnnnnnnnnnnnnnnnnnnnn 

□nnnnnnnnnnnnnnnnnnnnnnn 

□nnnnnnnnnnnnnnnnnnnnnnn 

□□□□□nnnnnnnnnnnnnnnnnnn 

□□□□□nnnnnnnnnnnnnnnnnnn 

□nnnnnnnnnnnnnnnnnnnnnnn 

□nnnnnnnnnnnnnnnnnnnnnnn 

□nnnnnnnnnnnnnnnnnnnnnnn 

□nnnnnnnnnnnnnnnnnnnnnnn 

□nnnnnnnnnnnnnnnnnnnnnnn 


1 

0.8 

0.6 

0.4 

0.2 

0 

- 0.2 

-0.4 

- 0.6 

- 0.8 


1 

0.8 

0.6 

0.4 

0.2 

0 

- 0.2 

-0.4 

- 0.6 

- 0.8 

r 


(d) BELPEX - Belgium 

for h,l = 0,..., 23 and selected markets from 17.12.2009 to 














































































































6 Figures 


15 


d-1,1 


d-1,1 






-o 

Oh 


nnnnnnnnnnnnnnnnnnnnnnnn 

n nn n nnnnnnnnnnnnnnnn n nn n 

□ □□□□□□□□□□□□□□□□□□□ nnnn 

nnnnnnnnnnnnnnnnnnnnnnnn 

n nn n nnnnnnnnnnnnnnnn n nn n 

nnnn nnnnnnnnnnnnnnnn nnnn 

n nn n nnnnnnnnnnnnnnnn n nn n 

□□□□□□□□□□□□□□□□□□□□□□□□ 

□□□□ □□□□nnnnnnnnnnnn nnnn 

nn n nnnnnnnnnnnnnnnn n nn n 

□□□□□□□□□□□□□□□□□□□□□□□ 

□□□□□□□□□□□nnnnnnnnnnnn 

□□□□□□□□□□nnnnnnnnnnnnn 

nn n nnnnnnnnnnnnnnnn n nn n 

nn n nnnnnnnnnnnnnnnn n nn n 

□□□□□□□□□□□□□□□□□□□□□□□ 

□□□□□□□□□□□□□□nnnnnnngg 

ggggggggggggggggnnnnnggg 

ggggggggggggggggggnnnggg 

□□□□ggggggggngggggnnnnnn 

n nn n nnnnnnnnnnnnnnnn n nn n 

□□□□ □□gggggggggggggg nnnn 

nnnnngggggggggggggggnnnn 


(a) APX - Netherlands 


1 

0.8 

0.6 

0.4 

0.2 

0 

- 0.2 

-0.4 

- 0.6 

- 0.8 

-r 




HtNrT!t'o^t'-ooo\o—'(Nro^ 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

•=11 
CV 12 

^ 13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 


nnnnnnnnaaannnnnnnnnnnn n 

nnnnnnnaannannnnnannnnnni 

nnnnnnnnnnannnnnnnnnnnnn 


nnnnnannnannnnnnnnnnnn 
□□□□□□□□□□□□□□□□□□□nnn 
nnnnnaanannnnnnnnnnnnn 
□□□□□nnnnnnnnnnnnnnnnnnn 

□□□□□nannngggggggggggg 
□□□□□□□□□□□ggggggggggg 
nnnnnnnnnnggggggggggga 
□nannnnnnaggggagggggaa 
□□□□□□□□□□□ggggggggggg 
□ggggggnnnaannnnnnnnnnnn 
ngggggnnnnnnnnnnnnnnnnnn 
□gggggnDDDDD nnnnnnnnnnnc 
□ggggggnnnnnnnnnnnnnnnnn 
nnnnnnnnnnnnnnnnnnnnnnnn 
□gggggnDDDDD □nnnnnnnnngc 
□□□□□□□DDDDPnnnnnnnnnnnn 
□□□□□□□DDDDPnnnnnnn nn nnn 
□□□□□□□DDDDPnnnnnnnn nnn n 
□□□□□□□DDDDPnnnnnnnnnnnn 
nnnnnnnDDDDDnnnnnnnnn nnn 
nnnnnnnDDDDDnnnnnnnnnnnn 


(b) Nordpool - Denmark West 
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(c) Nordpool - Denmark East 
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Figure 5: Sample correlations of Cor {Pd,h, Pd-i,i) for h, 
12.08.2014. 


I = 0,...,23 and selected markets from 17.12.2009 to 
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6.2 Weekly sample mean plots 
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Figure 6: Weekly sample mean of Pd^h for selected markets. 
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Figure 7: Weekly sample mean of Pd^h for selected markets. 
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6.3 


MAE/i plots 
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Figure 8: MAE/j for selected markets 
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6.4 


RMSEft plots 
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Figure 10: RMSE?i for selected markets 
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7 Tables of parameter importance 

EXAA Germany &; Austria 


Importance: 1 

Importance: 2 

Importance: 3 

Importance: 4 

Importance: 5 

Importance: 6 

Importance: 7 

0 

23@1 (65.3) 

3@I (10.2) 

Sun ( 6.6) 

19@I ( 4.2) 

Sat ( 3.7) 

401 ( 3.2) 

1801 ( 3.1) 

1 

23@1 (50.0) 

4@I (10.3) 

Sun ( 8.9) 

3@I ( 7.7) 

19@1 ( 6.3) 

2302 ( 2.5) 

1801 ( 2.0) 

2 

23@1 (39.1) 

3@I (13.2) 

Sun ( 8.1) 

23@2 ( 7.3) 

4@1 ( 6.9) 

1901 ( 4.9) 

2207 ( 3.1) 

3 

23@1 (31.9) 

3@I (10.6) 

4@I ( 9.7) 

23@2 ( 8.3) 

Sun ( 8.3) 

2207 ( 4.4) 

2001 ( 3.7) 

4 

23@1 (30.6) 

4@I (19.7) 

Sun (11.1) 

23@2 ( 5.8) 

2207 ( 5.5) 

406 ( 5.1) 

2001 ( 4.6) 

5 

23@1 (27.2) 

Sun (14.9) 

4@1 ( 7.5) 

20@I ( 6.7) 

507 ( 6.5) 

2207 ( 6.1) 

501 ( 4.6) 

6 

Sun (25.7) 

20@1 (17.6) 

23@1 (14.2) 

Sat (12.6) 

607 ( 7.8) 

Mon ( 5.8) 

2207 ( 4.5) 

7 

Sun (26.0) 

20@I (17.8) 

Sat (14.6) 

Mon ( 7.7) 

2301 ( 7.0) 

1801 ( 5.0) 

707 ( 4.8) 

8 

Sun (23.6) 

20@I (15.0) 

Sat (13.5) 

23@1 ( 8.8) 

Mon ( 8.1) 

701 ( 5.4) 

1801 ( 5.1) 

9 

Sun (23.8) 

20@I (15.4) 

Sat (13.5) 

Mon ( 9.8) 

2301 ( 8.0) 

1101 ( 5.1) 

1801 ( 5.0) 

10 

Sun (19.5) 

II@I (13.9) 

Sat (12.4) 

Mon (10.9) 

2001 ( 8.6) 

2301 ( 6.7) 

2201 ( 6.5) 

11 

Sun (14.2) 

Sat (11.5) 

Mon (11.3) 

11@1 (10.6) 

2201 ( 9.9) 

1201 ( 6.9) 

2301 ( 6.0) 

12 

12@1 (12.1) 

Mon (11.9) 

Sun (11.4) 

Sat ( 9.3) 

2301 ( 8.9) 

1301 ( 8.8) 

2201 ( 7.6) 

13 

13@1 (15.0) 

Sun (13.9) 

Mon (11.9) 

23@1 (10.1) 

Sat (10.0) 

1701 ( 6.2) 

1407 ( 5.9) 

14 

Sun (15.3) 

14@1 (15.0) 

Mon (11.2) 

23@1 (10.5) 

Sat ( 9.9) 

1701 ( 7.8) 

1407 ( 6.2) 

15 

Sun (17.3) 

15@1 (11.2) 

Mon (10.9) 

17@1 (10.6) 

Sat ( 9.4) 

2301 ( 8.7) 

1407 ( 5.9) 

16 

Sun (16.9) 

17@1 (15.2) 

Mon ( 8.6) 

Sat ( 8.3) 

2301 ( 7.7) 

1607 ( 6.8) 

1706 ( 6.2) 

17 

17@1 (32.2) 

Sun (13.0) 

Mon (11.1) 

17@6 (10.3) 

1707 ( 9.1) 

Sat ( 5.9) 

2207 ( 3.8) 

18 

18@1 (42.1) 

Mon (14.2) 

Sun ( 9.6) 

18@6 ( 9.4) 

1807 ( 8.8) 

Sat ( 4.5) 

1706 ( 2.8) 

19 

19@1 (35.2) 

Mon (11.3) 

Sun ( 8.8) 

20@1 ( 8.0) 

1906 ( 7.4) 

Sat ( 6.7) 

1806 ( 4.4) 

20 

20@1 (42.6) 

20@6 (11.5) 

Mon (10.0) 

Sat ( 9.1) 

Sun ( 6.5) 

2302 ( 2.7) 

2004 ( 2.3) 

21 

21@1 (38.5) 

21@6 (12.0) 

Mon ( 9.7) 

Sat ( 7.7) 

2201 ( 5.8) 

2302 ( 4.5) 

21035 ( 3.3) 

22 

22@1 (44.8) 

22@6 ( 9.2) 

23@2 ( 6.6) 

Sat ( 6.0) 

2301 ( 5.7) 

Mon ( 5.5) 

005 ( 5.4) 

23 

23@1 (36.8) 

22@1 ( 9.4) 

23@2 ( 7.5) 

23@6 ( 7.1) 

2304 ( 5.4) 

23028 ( 5.2) 

Mon ( 4.8) 


Table 4: Most relevant coefficients for each hour. H@L represents the estimated parameter for Yd-L,H for a model on 
Yd^h- Sun, Mon, ..Sat represent the weekday dummies. The corresponding th^i value is given in parenthesis. 



7 Tables of parameter importance 23 


EPEX Germany &; Austria 


Importance: 1 

Importance: 2 

Importance: 3 

Importance: 4 

Importance: 5 

Importance: 6 

Importance: 7 

0 

23@1 (70.6) 

19@I ( 5.9) 

I@I ( 5.7) 

18@1 ( 4.1) 

2101 ( 3.6) 

401 ( 3.1) 

Sat 

(2.2) 

1 

23@1 (67.3) 

19@I ( 8.6) 

3@I ( 7.1) 

4@1 ( 5.9) 

Sun ( 4.7) 

501 ( 2.6) 

601 

(2.0) 

2 

23@1 (56.3) 

4@I ( 9.4) 

19@I ( 8.9) 

3@1 ( 7.9) 

Sun ( 6.6) 

2207 ( 5.1) 

2025 

( 1-6) 

3 

23@1 (45.0) 

4@I (14.3) 

3@I (12.2) 

19@1 ( 8.5) 

Sun ( 7.3) 

2207 ( 4.3) 

306 

(2.5) 

4 

23@1 (44.8) 

4@I (13.8) 

19@I (10.3) 

Sun (10.2) 

301 ( 7.4) 

2207 ( 2.5) 

4025 

(2.4) 

5 

23@1 (34.1) 

Sun (15.5) 

5@I ( 9.7) 

4@1 ( 8.2) 

Sat ( 7.4) 

1901 ( 7.2) 

2207 

(3.8) 

6 

Sun (20.4) 

23@I (20.3) 

Sat (12.8) 

6@1 ( 9.4) 

Mon ( 8.7) 

2001 ( 6.0) 

2302 

(4.1) 

7 

Sun (23.9) 

Sat (15.9) 

23@I (I2.I) 

Mon ( 7.8) 

2001 ( 6.3) 

601 ( 5.3) 

1901 

(4.5) 

8 

Sun (28.3) 

Sat (15.9) 

20@I ( 8.5) 

23@1 ( 7.9) 

1801 ( 7.0) 

Mon ( 5.4) 

2101 

(5.2) 

9 

Sun (25.9) 

20@I (13.3) 

Sat (13.3) 

23@1 ( 7.0) 

2201 ( 6.3) 

Mon ( 6.1) 

1801 

(5.5) 

10 

Sun (21.8) 

Sat (11.9) 

22@I ( 9.6) 

20@1 ( 9.0) 

Mon ( 6.6) 

1301 ( 6.1) 

2301 

(5.6) 

11 

Sun (17.2) 

22@1 (12.3) 

Sat (11.7) 

13@1 ( 9.2) 

Mon ( 8.0) 

2301 ( 5.8) 

1107 

(4.1) 

12 

Sun (15.6) 

13@1 (12.0) 

Sat (10.7) 

22@1 (10.2) 

Mon ( 8.1) 

2301 ( 8.0) 

1701 

(4.9) 

13 

Sun (18.6) 

Sat (11.2) 

13@1 (11.1) 

22@1 (10.0) 

Mon ( 7.9) 

2301 ( 7.3) 

1701 

(6.8) 

14 

Sun (18.7) 

Sat (10.1) 

Mon ( 8.2) 

23@1 ( 8.2) 

1601 ( 8.1) 

2201 ( 7.6) 

1301 

(4.2) 

15 

Sun (19.8) 

16@1 (12.5) 

Sat ( 8.9) 

Mon ( 8.0) 

2301 ( 7.1) 

1706 ( 6.0) 

1701 

(4.9) 

16 

Sun (19.3) 

17@1 (12.5) 

Sat ( 8.8) 

16@1 ( 8.5) 

Mon ( 8.4) 

2301 ( 7.1) 

1706 

(6.7) 

17 

17@1 (26.4) 

Sun (16.1) 

Mon (10.1) 

17@6 (10.0) 

Sat ( 6.7) 

1707 ( 6.3) 

1801 

(6.0) 

18 

18@1 (33.8) 

Sun (14.1) 

Mon ( 9.9) 

18@6 ( 9.4) 

Sat ( 6.0) 

1802 ( 5.9) 

1805 

(4.6) 

19 

19@1 (34.0) 

Mon ( 9.4) 

19@6 ( 8.5) 

Sun ( 8.2) 

1902 ( 7.6) 

806 ( 5.0) 

1907 

(4.7) 

20 

20@1 (23.9) 

20@2 (10.4) 

Sat ( 9.5) 

Mon ( 7.8) 

2101 ( 7.5) 

2006 ( 7.5) 

Sun 

(7.0) 

21 

21@1 (22.8) 

Sat (10.9) 

21@2 ( 9.9) 

Mon ( 6.2) 

21028 ( 4.7) 

Sun ( 4.6) 

2301 

(4.6) 

22 

22@1 (30.6) 

22@2 ( 9.5) 

Sat ( 7.7) 

23@3 ( 6.5) 

22028 ( 6.1) 

2205 ( 5.1) 

22035 

(5.0) 

23 

23@1 (24.8) 

22@1 (17.5) 

23@3 (10.7) 

22@2 ( 9.6) 

Sat ( 5.7) 

005 ( 5.0) 

1207 

(4.4) 


Table 5: Most relevant coefficients for each hour. H@L represents the estimated parameter for Yd-L,H for a model on 
Yd^h- Sun, Mon, ..Sat represent the weekday dummies. The corresponding value is given in parenthesis. 
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EPEX Switzerland 


Importance: 1 

Importance: 2 

Importance: 3 

Importance: 4 

Importance: 5 

Importance: 6 

Importance: 7 

0 

23@1 (50.8) 

22@I (11.4) 

I@I ( 6.1) 

Sat ( 5.6) 

6@1 ( 5.2) 

5@1 ( 4.2) 

I8@l 

(3.7) 

1 

23@1 (47.6) 

5@1 (13.6) 

2@I (12.0) 

18@I ( 6.5) 

Sun ( 5.6) 

Sat ( 4.4) 

2@3 

(3.9) 

2 

23@1 (32.6) 

2@1 (12.7) 

4@I (10.6) 

Sun ( 8.4) 

18@1 ( 8.1) 

5@1 ( 5.4) 

3@7 

(3.5) 

3 

23@1 (30.7) 

4@1 (23.9) 

Sun ( 9.6) 

18@I ( 6.0) 

4@7 ( 4.3) 

4@4 ( 4.2) 

5@1 

(3.7) 

4 

23@1 (22.6) 

4@1 (20.8) 

Sun (12.0) 

5@I ( 6.7) 

22@1 ( 6.6) 

4@4 ( 5.1) 

19@5 

(4.8) 

5 

23@1 (25.9) 

Sun (19.2) 

5@I (13.4) 

Sat ( 6.9) 

5@7 ( 5.6) 

18@1 ( 4.9) 

4@1 

(4.0) 

6 

Sun (27.4) 

23@1 (20.8) 

Sat (12.0) 

20@I ( 9.5) 

22@1 ( 4.8) 

6@6 ( 4.8) 

5@1 

(4.7) 

7 

Sun (28.1) 

Sat (12.4) 

22@I (12.0) 

20@I ( 9.9) 

23@1 ( 8.5) 

18@1 ( 5.4) 

Mon 

(4.5) 

8 

Sun (22.8) 

22@1 (11.3) 

19@I (10.3) 

Sat ( 8.8) 

Mon ( 6.7) 

10@1 ( 6.5) 

8@6 

(5.6) 

9 

Sun (16.9) 

10@1 (10.7) 

22@I ( 9.8) 

19@I ( 8.2) 

Mon ( 7.7) 

Sat ( 7.0) 

9@2 

(6.5) 

10 

Sun (15.8) 

22@1 (14.2) 

Mon (lO.I) 

1I@I ( 9.7) 

19@1 ( 9.5) 

Sat ( 7.8) 

I0@1 

(6.2) 

11 

Sun ( 8.3) 

22@1 ( 6.9) 

10@2 ( 6.3) 

14@I ( 5.0) 

Sat ( 4.7) 

Mon ( 4.4) 

I3@7 

(4.3) 

12 

Sun ( 8.2) 

10@2 ( 7.6) 

22@I ( 7.2) 

13@7 ( 7.1) 

23@1 ( 5.1) 

13@1 ( 4.8) 

1201 

(4.5) 

13 

12@1 (16.8) 

Sun (15.8) 

22@I (11.5) 

13@7 ( 9.0) 

23@1 ( 9.0) 

Sat ( 8.8) 

Mon 

(8.4) 

14 

Sun (17.0) 

22@1 (14.5) 

Sat (10.2) 

Mon ( 9.9) 

13@7 ( 8.6) 

12@1 ( 8.3) 

I4@l 

(7.3) 

15 

Sun (12.7) 

22@1 (10.6) 

10@2 ( 6.9) 

Mon ( 6.7) 

Sat ( 6.2) 

23@1 ( 5.3) 

I3@7 

(5.3) 

16 

Sun (14.5) 

2201 (13.3) 

Mon ( 9.1) 

17@1 ( 8.9) 

Sat ( 6.8) 

14@1 ( 6.1) 

I8@l 

(5.1) 

17 

17@1 (15.0) 

18@1 (10.3) 

Sun ( 8.4) 

22@1 ( 6.5) 

10@2 ( 5.6) 

17@6 ( 4.7) 

Mon 

(4.7) 

18 

18@1 (40.6) 

18@6 ( 9.9) 

Mon ( 8.9) 

Sun ( 8.8) 

Sat ( 5.8) 

23@1 ( 4.8) 

22@1 

(4.4) 

19 

19@1 (28.9) 

10@2 (11.2) 

14@2 (11.0) 

Sun ( 5.2) 

22@1 ( 4.9) 

23@1 ( 4.7) 

I8@6 

(4.2) 

20 

20@1 (26.7) 

23@1 ( 7.7) 

Mon ( 6.8) 

22@1 ( 6.5) 

9@2 ( 5.9) 

15@2 ( 5.5) 

Sun 

(5.4) 

21 

22@1 (18.1) 

23@1 (17.7) 

21@1 (15.9) 

Mon ( 9.8) 

20@1 ( 9.7) 

13@7 ( 4.6) 

6@6 

(4.1) 

22 

22@1 (29.6) 

23@1 (22.1) 

Mon ( 9.0) 

13@6 ( 6.1) 

13@7 ( 4.7) 

Sat ( 4.0) 

20@1 

(3.8) 

23 

23@1 (39.6) 

22@1 (20.8) 

Mon ( 9.8) 

13@7 ( 6.3) 

6@6 ( 4.6) 

18@1 ( 4.6) 

0@2 

(4.3) 


Table 6: Most relevant coefficients for each hour. H@L represents the estimated parameter for Yd-L,H for a model on 
Yd^h- Sun, Mon, ..Sat represent the weekday dummies. The corresponding th^ value is given in parenthesis. 
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BELPEX Belgium 


Importance: 1 

Importance: 2 

Importance: 3 

Importance: 4 

Importance: 5 

Importance: 6 

Importance: 7 

0 

23@1 (29.3) 

0@I (10.4) 

2@I ( 8.8) 

7@7 ( 6.9) 

7@8 ( 5.1) 

101 ( 4.3) 

006 

(4.3) 

1 

2@1 (19.9) 

0@I (12.2) 

23@I ( 8.8) 

9@1 ( 6.4) 

1@1 ( 5.6) 

2101 ( 5.4) 

1301 

(5.3) 

2 

4@1 (23.2) 

2@I (22.9) 

0@I ( 7.6) 

I@1 ( 7.1) 

12@6 ( 6.2) 

1006 ( 5.9) 

1708 

(3.2) 

3 

4@1 (22.8) 

2@I (11.9) 

3@I ( 8.6) 

10@6 ( 5.4) 

0@1 ( 5.4) 

807 ( 4.2) 

101 

(3.6) 

4 

4@1 (33.7) 

5@1 (10.7) 

10@6 (10.4) 

3@I ( 8.5) 

0@1 ( 6.1) 

407 ( 5.0) 

101 

(4.2) 

5 

5@1 (38.7) 

4@1 (12.2) 

6@6 (11.5) 

5@7 ( 8.6) 

15@8 ( 5.1) 

1006 ( 5.1) 

602 

(4.1) 

6 

5@1 (40.0) 

6@7 (13.1) 

6@1 (11.2) 

7@6 ( 7.0) 

6@14 ( 5.9) 

606 ( 5.7) 

707 

(5.5) 

7 

5@1 (29.9) 

7@7 (16.1) 

7@1 (12.9) 

7@14 ( 9.9) 

7@6 ( 7.1) 

7021 ( 4.6) 

807 

(3.8) 

8 

7@1 (12.0) 

10@7 (11.5) 

5@1 (10.9) 

8@7 (10.6) 

8@14 ( 7.7) 

301 ( 6.5) 

201 

(6.2) 

9 

10@7 (59.7) 

4@1 (21.9) 

5@1 (18.3) 

0@1 ( 0.0) 

0@2 ( 0.0) 

003 ( 0.0) 

004 

(0.0) 

10 

10@7 (17.1) 

2@1 (13.5) 

10@14 ( 7.4) 

10@35 ( 6.2) 

4@1 ( 6.0) 

601 ( 5.5) 

lOOl 

(5.1) 

11 

2@1 (11.6) 

10@7 ( 9.1) 

11@14 ( 7.8) 

12@7 ( 6.6) 

11@35 ( 6.0) 

11021 ( 4.4) 

401 

(4.1) 

12 

2@1 (15.0) 

12@1 ( 8.9) 

12@14 ( 8.4) 

13@7 ( 7.4) 

12@35 ( 6.3) 

1007 ( 4.8) 

1207 

(3.7) 

13 

2@1 (20.3) 

13@7 (13.4) 

13@14 ( 8.6) 

3@1 ( 7.5) 

13@1 ( 7.0) 

1607 ( 6.8) 

13035 

(6.7) 

14 

3@1 (11.4) 

2@1 (10.6) 

16@7 ( 9.6) 

13@1 ( 8.1) 

14@14 ( 7.8) 

14021 ( 6.4) 

14035 

(5.3) 

15 

2@1 (11.3) 

15@14 (10.1) 

3@1 (10.0) 

16@7 ( 6.7) 

15@21 ( 5.9) 

1601 ( 5.6) 

15035 

(5.2) 

16 

2@1 (13.4) 

16@14 (10.6) 

17@1 ( 9.0) 

3@1 ( 8.1) 

16@7 ( 7.5) 

1707 ( 4.9) 

16028 

(4.7) 

17 

17@1 (14.2) 

17@7 (12.4) 

2@1 (10.4) 

17@14 ( 7.8) 

17@6 ( 5.7) 

401 ( 4.3) 

701 

(3.8) 

18 

18@1 (26.2) 

2@1 (11.0) 

18@7 ( 7.1) 

18@14 ( 6.5) 

17@7 ( 6.0) 

701 ( 5.6) 

2101 

(4.4) 

19 

19@1 (14.1) 

2@1 (13.4) 

19@7 ( 8.2) 

12@7 ( 7.1) 

7@1 ( 6.7) 

19014 ( 4.9) 

1902 

(4.5) 

20 

20@1 (19.0) 

2@1 (10.4) 

12@7 ( 5.6) 

20@2 ( 5.3) 

20@6 ( 5.3) 

301 ( 4.7) 

701 

(4.3) 

21 

2@1 (15.9) 

21@7 (15.3) 

21@1 (13.0) 

20@1 ( 7.4) 

21@6 ( 6.2) 

401 ( 5.7) 

21014 

(4.8) 

22 

2@1 (18.6) 

22@6 (11.4) 

22@1 ( 9.8) 

4@1 ( 8.9) 

2107 ( 7.9) 

2202 ( 4.8) 

lOOl 

(4.6) 

23 

2@1 (26.9) 

4@1 (17.3) 

23@1 ( 8.5) 

23@2 ( 7.0) 

2107 ( 5.5) 

1201 ( 5.4) 

1207 

(5.2) 


Table 7: Most relevant coefficients for each hour. H@L represents the estimated parameter for Yd-L,H for a model on 
Yd^h- Sun, Mon, ..Sat represent the weekday dummies. The corresponding value is given in parenthesis. 
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APX Netherlands 


Importance: 1 

Importance: 2 

Importance: 3 

Importance: 4 

Importance: 5 

Importance: 6 

Importance: 7 

0 

23@1 (52.2) 

2@I (14.7) 

Sat ( 5.8) 

9@I ( 5.3) 

2201 ( 5.0) 

2101 ( 2.7) 

108 

(2.1) 

1 

23@1 (22.9) 

2@I (12.0) 

5@I ( 7.2) 

16@I ( 5.3) 

106 ( 5.3) 

1901 ( 5.0) 

2201 

(4.7) 

2 

23@1 (22.5) 

2@I (20.8) 

5@I (11.1) 

19@I ( 7.8) 

301 ( 6.3) 

1601 ( 5.6) 

2201 

(5.3) 

3 

3@1 (20.3) 

23@I (16.5) 

5@1 ( 9.1) 

4@6 ( 5.6) 

Sun ( 4.6) 

401 ( 4.5) 

2101 

(4.2) 

4 

4@1 (12.3) 

3@I (12.3) 

23@1 (12.2) 

5@1 (II.5) 

Sun ( 9.3) 

507 ( 5.6) 

2101 

(5.6) 

5 

5@1 (23.5) 

Sun (14.9) 

23@1 (10.5) 

6@7 ( 9.9) 

2101 ( 8.1) 

1901 ( 6.9) 

506 

(6.7) 

6 

Sun (23.2) 

6@7 (12.5) 

Sat (12.4) 

23@I ( 7.0) 

1901 ( 6.9) 

501 ( 6.6) 

601 

(5.9) 

7 

Sun (27.5) 

Sat (17.3) 

20@1 ( 6.9) 

19@I ( 6.6) 

701 ( 6.1) 

Mon ( 5.5) 

707 

(3.8) 

8 

Sun (27.2) 

Sat (17.0) 

7@1 ( 6.9) 

7@7 ( 6.0) 

Mon ( 5.7) 

1901 ( 5.7) 

2001 

(5.7) 

9 

Sun (26.3) 

Sat (11.6) 

23@1 ( 5.2) 

20@I ( 4.3) 

1801 ( 4.1) 

707 ( 4.0) 

2101 

(3.7) 

10 

Sun (25.9) 

9@7 (10.8) 

12@1 (10.5) 

18@I ( 8.0) 

Sat ( 6.9) 

2201 ( 5.4) 

1006 

(4.8) 

11 

Sun (23.4) 

I2@I (11.8) 

Sat ( 9.6) 

22@I ( 7.0) 

Mon ( 5.4) 

1601 ( 5.2) 

1007 

(4.1) 

12 

Sun (14.3) 

12@1 (11.5) 

12@6 ( 7.6) 

Sat ( 6.1) 

Mon ( 5.8) 

2201 ( 4.2) 

12014 

(4.2) 

13 

Sun (16.2) 

12@1 (11.2) 

Mon (10.6) 

Sat (10.0) 

1601 ( 9.8) 

1206 ( 6.1) 

2201 

(5.6) 

14 

Sun (17.8) 

Sat (12.3) 

16@1 (10.1) 

1201 ( 8.9) 

Mon ( 8.5) 

1507 ( 7.3) 

2201 

(6.2) 

15 

Sun (18.7) 

16@1 (14.6) 

Sat (11.3) 

15@7 ( 9.0) 

Mon ( 6.5) 

2201 ( 6.5) 

1506 

(5.7) 

16 

Sun (16.7) 

16@1 (13.7) 

Sat (10.3) 

Mon ( 7.1) 

1607 ( 6.4) 

1702 ( 6.1) 

1701 

(5.5) 

17 

17@1 (27.3) 

Sun (11.5) 

17@7 (10.1) 

Sat ( 9.2) 

1706 ( 7.5) 

1702 ( 5.6) 

Mon 

(5.4) 

18 

18@1 (34.5) 

18@2 ( 8.2) 

Mon ( 7.9) 

18@7 ( 6.1) 

1803 ( 6.0) 

1706 ( 5.3) 

1806 

(4.4) 

19 

19@1 (44.1) 

19@2 (15.4) 

19@6 ( 8.6) 

19@3 ( 6.6) 

Mon ( 5.2) 

19014 ( 5.1) 

706 

(3.9) 

20 

20@1 (33.8) 

20@2 (16.3) 

20@6 ( 8.2) 

20@3 ( 7.8) 

20013 ( 5.2) 

1506 ( 5.1) 

606 

(4.3) 

21 

21@1 (25.4) 

21@2 (15.0) 

21@4 ( 6.0) 

21013 ( 5.9) 

Sat ( 5.3) 

2106 ( 5.1) 

OOl 

(4.7) 

22 

22@1 (21.3) 

22@2 (13.2) 

23@1 (12.7) 

22013 ( 5.4) 

2203 ( 4.7) 

Sat ( 4.4) 

2308 

(4.2) 

23 

23@1 (27.2) 

22@1 (15.0) 

23@2 ( 8.8) 

2202 ( 6.1) 

2303 ( 5.4) 

2308 ( 5.4) 

23028 

(5.0) 


Table 8: Most relevant coefficients for each hour. H@L represents the estimated parameter for Yd-L,H for a model on 
Yd^h- Sun, Mon, ..Sat represent the weekday dummies. The corresponding th^ value is given in parenthesis. 
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Nordpool Denmark West 


Importance: 1 

Importance: 2 

Importance: 3 

Importance: 4 

Importance: 5 

Importance: 6 

Importance: 7 

0 

23@1 (70.2) 

4@I ( 9.4) 

17@I ( 5.4) 

23@2 ( 4.4) 

4@2 ( 4.1) 

1501 ( 2.6) 

304 ( 2.4) 

1 

23@1 (61.8) 

4@1 (15.3) 

17@I ( 6.6) 

4@2 ( 4.3) 

3@4 ( 4.2) 

1606 ( 2.7) 

301 ( 1.8) 

2 

23@1 (62.0) 

4@I (I9.I) 

3@I ( 5.3) 

17@I ( 4.4) 

3@4 ( 2.7) 

Sun ( 2.4) 

406 ( 2.0) 

3 

23@1 (48.8) 

4@I (21.0) 

4@6 (11.1) 

17@I ( 3.9) 

3@1 ( 3.7) 

501 ( 3.7) 

1606 ( 3.2) 

4 

23@1 (48.3) 

4@1 (23.8) 

Sun ( 6.7) 

5@I ( 6.2) 

4@6 ( 4.7) 

301 ( 4.0) 

1701 ( 2.5) 

5 

23@1 (37.2) 

5@1 (16.1) 

Sun (12.8) 

4@I ( 9.9) 

Sat ( 6.9) 

607 ( 5.2) 

1801 ( 3.7) 

6 

Sun (16.0) 

23@1 (14.6) 

21@1 (13.9) 

6@1 (12.0) 

Sat (11.1) 

Mon ( 7.2) 

1901 ( 6.4) 

7 

13@1 (40.0) 

16@1 (11.2) 

18@1 (10.5) 

Sat (10.4) 

Sun ( 8.9) 

Mon ( 8.4) 

Fri ( 6.3) 

8 

13@1 (45.7) 

Sat (13.0) 

Sun (12.5) 

18@1 (11.1) 

Mon (10.2) 

Fri ( 7.4) 

OOl ( 0.0) 

9 

13@1 (41.2) 

16@1 (25.2) 

17@1 ( 8.2) 

Mon ( 6.7) 

Fri ( 5.1) 

1801 ( 4.7) 

Sat ( 4.0) 

10 

16@1 (33.1) 

13@1 (28.7) 

17@1 (13.0) 

15@1 (10.4) 

Mon ( 5.3) 

Fri ( 3.8) 

Sat ( 2.1) 

11 

16@1 (32.6) 

13@1 (29.8) 

17@1 (13.1) 

15@1 ( 9.6) 

Mon ( 5.7) 

Fri ( 3.5) 

Sat ( 2.4) 

12 

22@1 (13.0) 

21@1 (10.0) 

Sun ( 9.7) 

13@1 ( 8.8) 

Sat ( 8.1) 

1801 ( 7.8) 

Mon ( 6.0) 

13 

Sun (10.9) 

22@1 (10.4) 

21@1 (10.2) 

Sat ( 8.8) 

17@1 ( 7.0) 

1307 ( 6.3) 

Mon ( 6.0) 

14 

Sun (12.3) 

21@1 (10.4) 

Sat ( 9.1) 

17@1 ( 9.0) 

23@1 ( 7.7) 

14028 ( 7.0) 

2201 ( 6.5) 

15 

Sun (13.1) 

21@1 (10.1) 

17@1 ( 9.4) 

Sat ( 8.9) 

23@1 ( 8.6) 

15028 ( 6.7) 

1607 ( 6.6) 

16 

17@1 (16.9) 

16@7 (11.3) 

16@1 (10.8) 

23@1 (10.4) 

Sun ( 9.6) 

1706 ( 6.8) 

Mon ( 6.7) 

17 

17@1 (22.4) 

18@1 (12.1) 

16@7 ( 9.1) 

17@6 ( 8.8) 

Mon ( 8.7) 

Sun ( 8.2) 

Sat ( 5.9) 

18 

18@1 (40.1) 

Mon ( 8.0) 

18@2 ( 7.5) 

19@1 ( 7.3) 

17@6 ( 6.0) 

1806 ( 4.9) 

Sun ( 4.9) 

19 

19@1 (47.7) 

12@7 (10.0) 

19@2 ( 8.0) 

20@5 ( 5.6) 

18@6 ( 5.4) 

1406 ( 4.0) 

2001 ( 3.5) 

20 

21@1 (22.0) 

19@1 (10.8) 

20@1 ( 7.7) 

20@6 ( 7.6) 

19@2 ( 6.9) 

1207 ( 6.8) 

2005 ( 5.4) 

21 

21@1 (34.8) 

12@7 ( 7.0) 

23@1 ( 6.7) 

22@1 ( 6.6) 

21@5 ( 6.0) 

606 ( 5.7) 

2108 ( 4.7) 

22 

22@1 (31.6) 

22@5 (11.5) 

23@1 ( 8.2) 

22@2 ( 8.1) 

2101 ( 5.8) 

1207 ( 5.3) 

2204 ( 4.5) 

23 

23@1 (43.3) 

22@2 ( 8.8) 

23@3 ( 4.4) 

21@1 ( 4.4) 

2306 ( 3.5) 

504 ( 3.0) 

2205 ( 3.0) 


Table 9: Most relevant coefficients for each hour. H@L represents the estimated parameter for Yd-L,H for a model on 
Yd^h- Sun, Mon, ..Sat represent the weekday dummies. The corresponding value is given in parenthesis. 
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Nordpool Denmark East 


Importance: 1 

Importance: 2 

Importance: 3 

Importance: 4 

Importance: 5 

Importance: 6 

Importance: 7 

0 

23@1 (68.12) 

4@1 ( 8.53) 

4@2 ( 8.08) 

19@1 ( 4.19) 

0@1 ( 3.25) 

2302 ( 2.11) 

405 

( 1-73) 

1 

23@1 (64.05) 

4@1 (14.74) 

4@2 ( 7.74) 

3@1 ( 7.09) 

4@5 ( 3.61) 

1901 ( 2.44) 

406 

( 0.33) 

2 

23@1 (55.15) 

4@1 (17.31) 

3@1 (12.03) 

4@2 ( 5.01) 

4@5 ( 3.51) 

406 ( 3.10) 

Sun 

( 1.65) 

3 

23@1 (43.57) 

4@1 (22.45) 

3@1 (12.87) 

4@6 (11.09) 

4@2 ( 2.93) 

Sun ( 2.82) 

405 

( 2.73) 

4 

23@1 (42.37) 

4@1 (28.57) 

3@1 ( 8.57) 

4@6 ( 6.15) 

Sun ( 3.94) 

402 ( 3.01) 

405 

( 2.87) 

5 

23@1 (33.17) 

4@1 (15.31) 

Sun (11.31) 

5@1 (10.39) 

5@6 ( 5.68) 

Sat ( 5.64) 

1901 

( 5.58) 

6 

23@1 (19.80) 

Sun (13.82) 

6@1 (10.20) 

19@1 ( 9.52) 

Sat ( 9.12) 

2201 ( 6.43) 

2101 

( 6.40) 

7 

19@1 (23.52) 

23@1 (16.48) 

Sun (12.43) 

21@1 (11.61) 

Sat ( 9.63) 

1601 ( 8.89) 

Mon 

( 7.23) 

8 

19@1 (23.30) 

21@1 (17.47) 

23@1 (14.63) 

Sun (10.59) 

Mon ( 7.90) 

Sat ( 7.64) 

1601 

( 7.44) 

9 

19@1 (24.02) 

21@1 (23.68) 

23@1 (13.26) 

Mon (10.61) 

Sun ( 8.49) 

807 ( 5.58) 

Sat 

( 5.51) 

10 

21@1 (33.92) 

19@1 (18.90) 

23@1 (11.54) 

Mon (11.25) 

8@7 ( 6.11) 

1601 ( 5.55) 

2201 

( 5.41) 

11 

21@1 (22.86) 

16@1 (12.18) 

19@1 (10.54) 

23@1 ( 9.81) 

Mon ( 8.30) 

Sun ( 5.61) 

801 

( 4.87) 

12 

21@1 (20.38) 

16@1 (13.66) 

19@1 ( 8.47) 

Sun ( 6.78) 

Mon ( 6.72) 

Sat ( 6.44) 

2301 

( 6.30) 

13 

21@1 (23.12) 

14@1 ( 9.47) 

Mon ( 7.04) 

Sat ( 6.94) 

19@1 ( 6.93) 

Sun ( 6.88) 

1601 

( 6.79) 

14 

21@1 (21.38) 

16@1 (14.27) 

23@1 ( 9.69) 

Sun ( 7.40) 

Sat ( 6.98) 

1706 ( 6.70) 

801 

( 5.72) 

15 

21@1 (17.41) 

16@1 (15.43) 

23@1 ( 8.86) 

8@1 ( 8.04) 

Sun ( 6.61) 

Mon ( 5.94) 

Sat 

( 5.70) 

16 

16@1 (27.37) 

23@1 ( 8.70) 

21@1 ( 7.92) 

Mon ( 6.52) 

1706 ( 6.11) 

1707 ( 5.18) 

1606 

( 4.64) 

17 

16@6 (34.01) 

17@1 (24.75) 

16@1 (24.74) 

16@7 (13.44) 

2301 ( 2.41) 

Sun ( 0.64) 

OOl 

( 0.00) 

18 

16@6 (40.52) 

19@1 (30.44) 

16@1 (19.99) 

16@7 ( 9.05) 

OOl ( 0.00) 

002 ( 0.00) 

003 

( 0.00) 

19 

19@1 (39.15) 

22@1 ( 9.71) 

23@1 ( 6.63) 

19@2 ( 6.60) 

Mon ( 6.52) 

1706 ( 5.66) 

1607 

( 3.73) 

20 

20@1 (24.92) 

22@1 (14.63) 

21@5 ( 6.47) 

19@1 ( 6.00) 

2301 ( 5.52) 

Sat ( 5.08) 

1706 

( 4.95) 

21 

21@1 (22.91) 

22@1 (15.72) 

23@1 (11.11) 

21@5 ( 9.02) 

1707 ( 5.04) 

Sat ( 4.53) 

Fri 

( 3.43) 

22 

22@1 (37.78) 

23@1 (17.21) 

22@2 (11.15) 

22@5 (10.50) 

1406 ( 3.15) 

2206 ( 2.97) 

2207 

( 2.77) 

23 

23@1 (52.37) 

22@2 ( 6.60) 

0@1 ( 5.79) 

22@5 ( 5.78) 

2201 ( 5.64) 

005 ( 3.39) 

303 

( 3.17) 


Table 10: Most relevant coefficients for each hour. H@L represents the estimated parameter for Yd-L,H for a model 
on Yd^h- Sun, Mon, ..., Sat represent the weekday dummies. The corresponding value is given in parenthesis. 
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Nordpool Sweden(4) 


Importance: 1 

Importance: 2 

Importance: 3 

Importance: 4 

Importance: 5 

Importance: 6 

Importance: 7 

0 

23@1 

(74.57) 

1@1 (10.65) 

4@1 ( 9.78) 

1803 ( 2.33) 

1702 ( 0.96) 

704 ( 0.63) 

Sat 

( 0.50) 

1 

23@1 

(46.44) 

4@1 (13.77) 

1@1 ( 9.20) 

2206 ( 5.28) 

2202 ( 4.92) 

303 ( 3.10) 

2102 

( 2.90) 

2 

23@1 

(41.80) 

4@1 (16.20) 

22@2 ( 6.91) 

2206 ( 6.81) 

301 ( 6.64) 

303 ( 4.90) 

2203 

( 3.90) 

3 

23@1 

(35.98) 

4@1 (22.09) 

22@6 ( 8.95) 

2202 ( 6.30) 

2203 ( 4.94) 

303 ( 4.83) 

405 

( 3.39) 

4 

23@1 

(27.78) 

4@1 (16.84) 

22@2 ( 5.81) 

2306 ( 5.18) 

405 ( 4.80) 

2206 ( 3.42) 

506 

( 2.89) 

5 

23@1 

(26.23) 

4@1 (17.77) 

2101 (10.58) 

2206 ( 8.55) 

2207 ( 5.45) 

507 ( 4.66) 

506 

( 4.40) 

6 

21@1 

(25.42) 

23@1 (14.07) 

2206 ( 8.11) 

807 ( 7.58) 

1107 ( 6.84) 

1101 ( 6.48) 

1901 

( 5.70) 

7 

21@1 

(17.21) 

8@7 (16.37) 

1107 ( 9.53) 

1101 ( 9.26) 

1007 ( 5.11) 

1601 ( 5.00) 

1801 

( 4.56) 

8 

21@1 

(24.41) 

8@7 (15.96) 

1107 (12.62) 

1101 ( 7.48) 

1901 ( 4.13) 

2103 ( 3.76) 

1601 

( 3.73) 

9 

21@1 

(30.62) 

8@7 (16.67) 

1107 (12.07) 

1901 ( 4.19) 

1101 ( 4.03) 

2103 ( 4.00) 

2007 

( 2.73) 

10 

21@1 

(26.27) 

8@7 (17.43) 

1107 (10.60) 

1801 ( 6.24) 

1101 ( 5.22) 

1007 ( 3.11) 

701 

( 3.06) 

11 

21@1 

(30.07) 

8@7 (13.90) 

1107 ( 9.53) 

1101 ( 5.41) 

607 ( 4.19) 

1801 ( 3.62) 

1601 

( 3.40) 

12 

21@1 

(19.11) 

22@1 (11.14) 

1601 ( 7.47) 

807 ( 5.79) 

1101 ( 5.62) 

1901 ( 4.79) 

702 

( 4.32) 

13 

21@1 

(17.35) 

8@7 ( 7.06) 

2201 ( 5.96) 

702 ( 5.10) 

801 ( 4.23) 

2301 ( 3.94) 

1801 

( 3.45) 

14 

21@1 

(18.97) 

8@7 ( 5.70) 

702 ( 4.75) 

1801 ( 4.40) 

1601 ( 4.23) 

801 ( 4.14) 

1101 

( 3.88) 

15 

21@1 

(20.49) 

16@1 ( 8.22) 

801 ( 7.50) 

807 ( 4.90) 

1107 ( 4.86) 

2301 ( 4.70) 

2201 

( 4.44) 

16 

21@1 

(13.72) 

16@1 (12.11) 

1101 ( 5.77) 

702 ( 5.64) 

1107 ( 4.82) 

2206 ( 4.33) 

2201 

( 4.25) 

17 

21@1 

(18.05) 

8@7 (12.49) 

1701 (10.94) 

901 (10.56) 

1107 (10.54) 

1101 ( 4.33) 

2103 

( 3.53) 

18 

21@1 

(19.43) 

8@7 (17.33) 

1107 (10.77) 

1901 ( 8.85) 

1801 ( 6.29) 

1101 ( 5.24) 

2103 

( 4.49) 

19 

19@1 

(11.97) 

21@1 (10.90) 

807 ( 9.15) 

1107 ( 8.18) 

1101 ( 7.31) 

2001 ( 5.64) 

705 

( 5.00) 

20 

21@1 

(23.22) 

8@7 ( 9.37) 

1107 ( 8.54) 

2001 ( 7.37) 

2201 ( 5.51) 

1901 ( 5.50) 

705 

( 4.95) 

21 

21@1 

(54.22) 

22@1 (26.19) 

2301 ( 7.75) 

1406 ( 3.83) 

Mon ( 2.46) 

1606 ( 1.92) 

Fri 

( 1-51) 

22 

22@1 

(66.34) 

23@1 (10.96) 

OOl ( 4.78) 

2205 ( 3.76) 

2101 ( 2.87) 

904 ( 2.06) 

2106 

( 1.94) 

23 

23@1 

(54.20) 

22@1 (26.52) 

2302 ( 6.36) 

OOl ( 3.83) 

103 ( 2.82) 

405 ( 2.70) 

2305 

( 1-72) 


Table 11: Most relevant coefficients for each hour. H@L represents the estimated parameter for for a model 

on Yd^h- Sun, Mon, ..., Sat represent the weekday dummies. The corresponding value is given in parenthesis. 
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POLPX Poland 


Importance: 1 

Importance: 2 

Importance: 3 

Importance: 4 

Importance: 5 

Importance: 6 

Importance: 7 

0 

23@1 (23.8) 

0@1 (13.9) 

7@1 ( 8.5) 

1@1 ( 7.0) 

2302 ( 6.8) 

007 ( 6.0) 

2201 

(3.5) 

1 

23@1 (28.1) 

1@1 (23.1) 

1@7 ( 9.6) 

3@1 ( 8.9) 

2201 ( 5.8) 

Sat ( 5.6) 

601 

(4.5) 

2 

23@1 (23.6) 

3@1 (15.2) 

2@1 (10.4) 

22@1 ( 8.9) 

307 ( 8.1) 

401 ( 7.7) 

501 

(6.7) 

3 

4@1 (23.3) 

23@1 (17.5) 

3@1 (15.3) 

3@7 ( 9.5) 

2201 ( 7.7) 

Sun ( 7.2) 

Sat 

(4.2) 

4 

4@1 (29.4) 

23@1 (13.2) 

5@1 (11.6) 

Sun (10.4) 

307 ( 8.6) 

2201 ( 8.2) 

405 

(4.6) 

5 

5@1 (32.6) 

Sun (10.9) 

22@1 (10.6) 

23@1 (10.1) 

507 ( 8.7) 

605 ( 5.1) 

2202 

(3.4) 

6 

6@1 (23.2) 

Sun (11.4) 

23@1 ( 9.0) 

6@7 ( 8.2) 

2201 ( 7.7) 

Mon ( 4.8) 

606 

(3.8) 

7 

Sun (17.8) 

23@1 (14.5) 

7@1 (14.2) 

7@7 (12.0) 

Mon ( 7.4) 

Sat ( 6.7) 

706 

(5.3) 

8 

8@1 (17.9) 

Sun (12.9) 

22@1 (11.3) 

8@7 (10.5) 

Mon ( 5.7) 

2208 ( 4.1) 

Sat 

(4.0) 

9 

9@1 (20.1) 

Sun (16.6) 

22@1 ( 7.5) 

Mon ( 7.2) 

Sat ( 5.3) 

907 ( 4.8) 

9014 

(4.5) 

10 

10@1 (20.4) 

Sun (17.1) 

Mon (10.1) 

21@1 ( 6.8) 

2201 ( 6.1) 

1201 ( 5.2) 

806 

(3.9) 

11 

12@1 (19.7) 

14@7 ( 9.7) 

Sun ( 7.4) 

14@8 ( 5.6) 

2201 ( 4.7) 

1203 ( 4.0) 

Mon 

(3.9) 

12 

12@1 (24.1) 

14@7 ( 9.1) 

12@21 ( 5.9) 

12014 ( 5.2) 

1308 ( 4.7) 

Mon ( 4.1) 

Sun 

(4.0) 

13 

12@1 (19.1) 

14@7 ( 8.3) 

13@14 ( 5.6) 

13021 ( 5.3) 

1502 ( 4.0) 

Mon ( 3.9) 

1308 

(3.8) 

14 

14@1 (12.3) 

14@7 ( 9.8) 

12@1 ( 9.7) 

1108 ( 6.5) 

14021 ( 6.5) 

14014 ( 5.1) 

1501 

(4.7) 

15 

15@1 (26.2) 

14@7 (11.5) 

Sun ( 6.7) 

1408 ( 6.3) 

15014 ( 5.5) 

Mon ( 5.4) 

2201 

(5.3) 

16 

16@1 (14.5) 

16@7 ( 8.9) 

17@1 ( 7.8) 

2201 ( 7.2) 

Sun ( 6.5) 

16021 ( 4.9) 

1501 

(4.7) 

17 

17@1 (41.4) 

17@6 (10.2) 

Sun ( 9.0) 

1801 ( 7.9) 

1607 ( 7.7) 

Mon ( 7.5) 

1608 

(5.6) 

18 

18@1 (51.5) 

18@6 ( 8.6) 

Sun ( 7.0) 

Mon ( 6.5) 

1607 ( 4.8) 

1608 ( 4.0) 

1808 

(3.8) 

19 

19@1 (46.2) 

18@8 ( 7.1) 

18@6 ( 5.1) 

Sun ( 4.2) 

Mon ( 4.1) 

1907 ( 3.6) 

2001 

(3.4) 

20 

20@1 (56.0) 

Mon ( 5.5) 

20@6 ( 4.6) 

2208 ( 4.5) 

2207 ( 3.9) 

Sun ( 2.9) 

20020 

(2.6) 

21 

21@1 (43.1) 

22@7 ( 8.6) 

20@1 ( 5.6) 

1108 ( 5.4) 

2106 ( 5.2) 

2107 ( 5.2) 

Sat 

(4.8) 

22 

22@1 (34.6) 

22@7 (20.7) 

22@8 ( 6.7) 

2206 ( 5.5) 

22035 ( 4.5) 

22021 ( 4.1) 

2301 

(3.8) 

23 

23@1 (51.4) 

23@7 (10.8) 

23@6 ( 9.1) 

Sat ( 5.6) 

Mon ( 5.2) 

23021 ( 4.2) 

103 

(2.4) 


Table 12: Most relevant coefficients for each hour. H@L represents the estimated parameter for Yd-L,H for a model 
on Yd^h- Sun, Mon, ..., Sat represent the weekday dummies. The corresponding value is given in parenthesis. 
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OTE Czech 


Importance: 1 

Importance: 2 

Importance: 3 

Importance: 4 

Importance: 5 

Importance: 6 

Importance: 7 

0 

23@1 (55.3) 

21@1 (10.1) 

4@1 ( 9.2) 

7@1 ( 8.2) 

Sat ( 5.1) 

Sun ( 3.9) 

3@1 

(3.1) 

1 

23@1 (43.4) 

4@1 (14.1) 

Sun (11.7) 

21@1 ( 9.7) 

7@1 ( 9.4) 

3@1 ( 4.5) 

19@1 

(2.4) 

2 

23@1 (37.3) 

4@1 (18.5) 

Sun (10.4) 

20@1 ( 7.3) 

18@1 ( 4.7) 

3@6 ( 4.5) 

19@1 

(4.4) 

3 

23@1 (33.3) 

4@1 (22.4) 

Sun (11.9) 

20@1 ( 9.2) 

19@1 ( 4.5) 

3@1 ( 4.1) 

18@1 

(4.1) 

4 

23@1 (28.4) 

4@1 (21.6) 

Sun (15.3) 

20@1 ( 8.3) 

3@6 ( 5.1) 

18@1 ( 4.5) 

19@1 

(4.2) 

5 

23@1 (25.5) 

Sun (20.2) 

20@1 (11.1) 

4@1 ( 9.4) 

5@6 ( 7.9) 

19@1 ( 4.2) 

6@7 

(3.1) 

6 

Sun (28.6) 

23@1 (16.8) 

20@1 (15.4) 

Sat (12.5) 

Mon ( 7.6) 

6@1 ( 5.4) 

19@1 

(4.5) 

7 

Sun (25.7) 

20@1 (17.3) 

Sat (14.1) 

23@1 ( 8.1) 

Mon ( 7.0) 

7@1 ( 4.6) 

18@1 

(4.4) 

8 

Sun (25.4) 

20@1 (16.4) 

Sat (13.7) 

Mon ( 7.7) 

8@1 ( 6.7) 

23@1 ( 5.9) 

8@7 

(3.9) 

9 

Sun (22.4) 

20@1 (18.4) 

Sat (11.0) 

11@1 ( 9.8) 

Mon ( 9.5) 

17@1 ( 6.9) 

8@7 

(6.4) 

10 

Sun (17.5) 

20@1 (12.0) 

11@1 (10.5) 

Mon (10.4) 

22@1 ( 9.4) 

Sat ( 9.1) 

11@7 

(5.3) 

11 

22@1 (14.1) 

Sun (12.6) 

Mon (11.3) 

11@1 (11.2) 

Sat ( 8.3) 

20@1 ( 6.3) 

12®! 

(5.7) 

12 

Sun (13.0) 

22@1 (11.1) 

Mon (10.3) 

14@1 ( 9.7) 

Sat ( 8.0) 

12@1 ( 6.6) 

14@7 

(4.5) 

13 

14@1 (14.3) 

Sun (14.1) 

22@1 (13.3) 

Mon (12.1) 

Sat ( 8.6) 

14@7 ( 7.1) 

17@1 

(5.4) 

14 

Sun (16.3) 

14@1 (12.2) 

Mon (10.6) 

22@1 ( 9.7) 

Sat ( 8.2) 

14@7 ( 7.2) 

16@1 

(6.9) 

15 

Sun (17.7) 

16@1 (14.5) 

Mon (11.1) 

22@1 ( 9.9) 

Sat ( 7.6) 

14@7 ( 6.8) 

17@1 

(6.7) 

16 

Sun (16.6) 

16@1 (14.9) 

17@1 (12.3) 

Mon (10.3) 

22@1 ( 9.4) 

16@7 ( 7.1) 

Sat 

(6.7) 

17 

17@1 (35.0) 

Sun (12.7) 

Mon (12.0) 

17@7 ( 9.5) 

17@6 ( 5.7) 

Sat ( 5.3) 

22@1 

(3.8) 

18 

18@1 (30.6) 

Mon (10.9) 

18@6 (10.5) 

Sun ( 9.4) 

18@7 ( 7.7) 

17@1 ( 7.1) 

Sat 

(4.9) 

19 

19@1 (39.4) 

Mon (12.0) 

19@6 ( 7.5) 

19@7 ( 6.9) 

20@1 ( 6.7) 

Sun ( 6.1) 

Sat 

(5.8) 

20 

20@1 (45.8) 

Mon (13.2) 

20@6 (10.9) 

Sat ( 6.5) 

20@7 ( 6.3) 

22@1 ( 4.2) 

20@21 

(3.4) 

21 

21@1 (30.6) 

22@1 (14.0) 

21@6 (10.1) 

Mon ( 9.5) 

Sat ( 6.2) 

21@7 ( 5.6) 

20@1 

(4.9) 

22 

22@1 (45.7) 

22@6 ( 8.6) 

Mon ( 5.6) 

2204 ( 5.5) 

Sat ( 4.2) 

22@21 ( 3.9) 

22@7 

(3.8) 

23 

22@1 (30.6) 

23@1 (26.1) 

23@6 ( 7.9) 

Mon ( 6.9) 

22@2 ( 4.1) 

23@4 ( 3.7) 

23@35 

(2.7) 


Table 13: Most relevant coefficients for each hour. H@L represents the estimated parameter for Yd-L,H for a model 
on Yd^h- Sun, Mon, ..., Sat represent the weekday dummies. The corresponding value is given in parenthesis. 



